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The lectures review the state of affairs in modern branch of mathematical 
physics called probabilistic topology. In particular we consider the following 
problems: (i) We estimate the probability of a trivial knot formation on the 
lattice using the Kauffman algebraic invariants and show the connection of this 
problem with the thermodynamic properties of 2D disordered Potts model; (ii) 
We investigate the limit behavior of random walks in multi-connected spaces 
and on non-commutative groups related to the knot theory. We discuss the 
application of the above mentioned problems in statistical physics of polymer 
chains. On the basis of non-commutative probability theory we derive some 
new results in statistical physics of entangled polymer chains which unite 
rigorous mathematical facts with more intuitive physical arguments. 
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I. INTRODUCTION 

It wouldn't be an exaggeration to say that contemporary physical science is becoming 
more and more mathematical. This fact is too strongly manifested to be completely ignored. 
Hence I would permit myself to bring forward two possible conjectures: 

(a) On the one hand there are hardly discovered any newly physical problem which would 
be beyond the well established methods of the modern theoretical physics. This leads to 
the fact that nowadays real physical problems seem to be less numerous than mathematical 
methods of their investigation. 

(b) On the other hand the mathematical physics is a fascinating field which absorbs new 
ideas from different branches of modern mathematics, translates them into the physical lan- 
guage and hence fills the abstract mathematical constructions by the new fresh content. This 
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ultimately leads to creating new concepts and stimulates seeking for newel deep conformities 
to natural laws in known physical phenomena. 

The penetration of new mathematical ideas in physics has sometimes rather paradoxical 
character. It is not a secret that difference in means (in languages) and goals of physicists 
and mathematicians leads to mutual misunderstanding, making the very subject of investi- 
gation obscure. What is true for general is certainly true for particular. To clarify the point, 
let us turn to statistics of entangled uncrossible random walks — the well-known subject of 
statistical physics of polymers. Actually, since 1970s, after Conway's works, when the first 
algebraic topological invariants — Alexander polynomials — became very popular in mathe- 
matical literature, physicists working in statistical topology have acquired a much more 
powerful topological invariant than the simple Gauss linking number. The constructive uti- 
lization of algebraic invariants in statistical physics of macromolecules has been developed 
in the classical works of A. Vologodskii and M. Frank-Kamenetskii However until re- 
cently in overwhelming majority of works the authors continue using the commutative Gauss 
invariant invariant just making references to its imperfectness. 

One of the reasons of such inertia consists in the fact that new mathematical ideas are 
often formulated as "theorems of existence" and it takes much time to retranslate them into 
physically acceptable form which may serve as a real computational tool. 

We intend to use some recent advances in algebraic topology and theory of random walks 
on non-commutative groups for reconsidering the old problem — evaluating of the entropy of 
randomly generated knots and entangled random walks in a given homotopic state. Let us 
emphasize that this is a real physical paper, and when it is possible the rigorous statements 
are replaced by some physically justified conjectures. Generally speaking, the work is devoted 
to an analysis of probabilistic problems in topology and their applications in statistical 
physics of polymer systems with topological constraints. 

Let us formulate briefly the main results of our work. 

1. The probability for a long random walk to form randomly a knot with specific topo- 
logical invariant is computed. This problem is considered using the Kauffman algebraic 
invariants and the connection with the thermodynamic properties of 2D Potts model with 
"quenched" and "annealed" disorder in interaction constants is discussed. 

2. The limit behavior of random walks on the non-commutative groups related to the 
knot theory is investigated. Namely, the connection between the limit distribution for the 
Lyapunov exponent of products of non-commutative random matrices — generators of "braid 
group" — and the asymptotic of powers ("knot complexity") of algebraic knot invariants is 
established. This relation is applied for calculating the knot entropy. In particular, it 
is shown that the "knot complexity" corresponds to the well known topological invariant, 
"primitive path" , repeatedly used is statistics of entangled polymer chains. 

3. The random walks on multi-connected manifolds is investigated using conformal 
methods and the nonabelian topological invariants are constructed. It is shown that many 
nontrivial properties of limit behavior of random walks with topological constraints can be 
explained in context of random walks on hyperbolic groups. 

Usage of the limit behavior of entangled random paths established above for investi- 
gation of the statistical properties of so-called "crumpled globule" (trivial ring without 
self- intersections in strongly contracted state). 

The connection between all these problems is shown in Table 1. 



3 



Spectral and quantum 
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FIG. 1. Links between topologically-probabilistic problems. Solid boxes - problems, discussed 
in the paper; dashed boxes - problems not included in the consideration. 



II. KNOT DIAGRAMS AS DISORDERED SPIN SYSTEMS 
A. Brief review of statistical problems in topology 

The interdependence of such branches of modern theoretical and mathematical physics 
as theory of integrable systems, algebraic topology and conformal field theory has proved to 
be a powerful catalyst of development of the new direction in topology, namely, of analytical 
topological invariants construction by means of exactly solvable statistical models. 

Today it is widely believed that the following three cornerstone findings have brought 
the fresh stream in topology: 

— It has been found the deep relation between the Temperley-Lieb algebra and the Hecke 
algebra representation of the braid group. This fact resulted in the remarkable geometrical 
analogy between the Yang-Baxter equations, appearing as necessary condition of the transfer 
matrix commutativity in the theory of integrable systems on the one hand, and one of 
Reidemeister moves, used in the knot invariant construction on the other hand. 

— It has been discovered that the partition function of the Wilson loop with the Chern- 
Simons action in the topological field theory coincides with the representation of the known 
nonabelian algebraic knot invariants written in terms of the time-ordered path integral. 

— The need for new solutions of the Yang-Baxter equations has given a power impetus 
to the theory of quantum groups. Later on the related set of problems was separated in the 
independent branch of mathematical physics. 

Of course the above mentioned findings do not exhaust the list of all brilliant achieve- 
ments in that field during the last decade, but apparently these new accomplishments have 
used profound "ideological" changes in the topological science: now we can hardly consider 
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topology as an independent branch of pure mathematics where each small step forward takes 
so much effort that it seems incidental. 

Thus in the middle of the 80s the "quantum group" gin was released. It linked by 
common mathematical formalism classical problems in topology, statistical physics and field 
theory. A new look at the old problems and the beauty of the formulated ideas made an 
impression on physicists and mathematicians. As a result, in a few last years the number 
of works devoted to the search of the new applications of the quantum group apparatus is 
growing exponentially going beyond the framework of original domains. As an example of 
persistent penetrating of the quantum group ideas in physics we can name the works on anyon 
superconductivity 0, intensively discussing problems on "quantum random walks" 0, the 
investigation of spectral properties of "quantum deformations" of harmonic oscillators @] 
and so on. 

The time will show whether such "quantum group expansion" is physically justified or 
it merely does tribute to today's fashion. However it is clear that physics has acquired new 
convenient language allowing to construct new "nonabelian objects" and to work with them. 

Among the vast amount of works devoted to different aspects of the theory of integrable 
systems, their topological applications connected to the construction of knot and link in- 
variants and their representation in terms of partition functions of some known 2D-models 
deserve our special attention. There exist several reviews ^ and books on that subject 
and our aim by no means consists in re-interpretation or compilation of their contents. We 
make an attempt of consecutive account of recently solved probabilistic problems in topology 
as well as attract attention to some interesting, still unsolved, questions lying on the border 
of topology and the probability theory. Of course we employ the knowledges acquired in the 
algebraic topology utilizing the construction of new topological invariants done by V.F.R. 
Jones 1^ and L.H. Kauffman 

Besides the traditional fundamental topological issues concerning the construction of new 
topological invariants, investigation of homotopic classes and fibre bundles we mark a set of 
adjoint but much less studied problems. First of all, we mean the problem of so-called "knot 
entropy" calculation. Most generally it can be formulated as follows. Take the lattice 
embedded in the space R^. Let be the ensemble of all possible closed nonselfintersecting 
A^-step loops with one common fixed point onZ^; by cu we denote the particular trajectory 
configuration. The question is: what is the probability Vn of the fact that the trajectory 
uj G Qn belongs to some specific homotopic class. Formally this quantity can be represented 
in the following way 

VN{lnv} = 7^ E ^ [Inv{a;} - Inv] = 

1 (II.l) 

E A [Inv{ri, . . . , rjv} - Inv] (l - A [r, - r,]) A [r/ ^ 



{ri,...,riv} 



-N 



where Inv{a;} is the functional representation of the knot invariant corresponding to the 
trajectory with the bond coordinates {ri . . . ,r7v}; Inv is the topological invariant character- 
izing the knot of specific homotopic type and A(a;) is the Kronecker function: A(x = 0) = 1 
and A(a; 7^ 0) = 0. The first A-function in Eq.( [ll.lD cuts the set of trajectories with the 



fixed topological invariant while the second and the third A-functions ensure the A^-step 
trajectory to be nonselfintersecting and to form a closed loop respectively. 
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The distribution function P7v{Inv} satisfies the normahzation condition 

J2 Pjv{Inv} = 1 (II.2) 

all liomotopic classes 

The entropy S'Arjlnv} of the given homotopic state of the knot represented by A^-step 
closed loop on reads 

5jv{Inv} = In [^]Jv^7v{Inv}] (11.3) 

The problem concerning the knot entropy determination has been discussed time and 
again by the leading physicists. However the number of new analytic results in that field was 
insufficient till the beginning of the 80s: in about 90 percents of published materials their 
authors used the Gauss linking number or some of its abelian modifications for classification 
of a topological state of knots and links while the disadvantages of this approach were 
explained in the rest 10 percent of the works. We do not include in this list the celebrated 
investigations of A. V. Vologodskii et al |I| devoted to the first fruitful usage of the nonabelian 
Alexander algebraic invariants for the computer simulations in the statistical biophysics. We 
discuss physical applications of these topological problems at length in Section 5. 

Despite of the clarity of geometrical image, the topological ideas are very hard to formal- 
ize because of the non-local character of topological constraints. Besides, the main difficulty 
in attempts to calculate analytically the knot entropy is due to the absence of convenient 
analytic representation of the complete topological invariant. Thus, to succeed, at least 
partially, in the knot entropy computation we simplify the general problem replacing it by 
the problem of calculating the distribution function for the knots with defined topological 
invariants. That problem differs from the original one because none of the known topo- 
logical invariants (Gauss linking number, Alexander, Jones, HOMFLY) are complete. The 
only exception is Vassiliev invariants |]^, which are beyond the scope of the present book. 
Strictly speaking we are unable to estimate exactly the correctness of such replacement of 
the homotopic class by the mentioned topological invariants. Thus under the definition of 
the topological state of the knot or entanglement we simply understand the determination 
of the corresponding topological invariant. 

The problems where uj (see Eq.( p.l.l| )) is the set of realizations of the random walk, i.e. 
the Markov chain are of special interest. In that case the probability to find a closed A^-step 
random walk in in some prescribed topological state can be presented in the following 
way 




...,r^}-Inv]5[r^] (11.4) 



where g (rj+i — Vj) is the probability to find j + 1th step of the trajectory in the point rj+i 
if jth step is in r^. In the limit a ^ and N oo {Na = L = const) in three-dimensional 
space we have the following expression for g (r^+i — r^) 

f W 3 / 3(r,^i-r,)^ \ / 3 fs ( dris)]') 
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Introducing the "time", s, along the trajectory we rewrite the distribution function PAr{Inv} 
(Eq. (|II.4| )) in the path integral form with the Wiener measure density 

VNilnv} = ^J---J ^{^} £ 5[Inv{r(s)} - Inv] (11.6) 



and the normalization condition is as follows 



Z= J2 ^7v{Inv} 



all different 
knot invariants 



The form of Eq. ( [11. 61 ) up to the Wick turn and the constants coincides with the scattering 



amplitude a of a free quantum particle in the multi-connected phase space. Actually, for 
the amplitude a we have 

a~ ^ exp I ^ f r\s)ds] (11.7) 



h 



all paths from given 
topological class 



If phase trajectories can be mutually transformed by means of continuous deformations, then 
the summation in Eq. ( |11.71) should be extended to all available paths in the system, but if 
the phase space consists of different topological domains, then the summation in Eq.( [11.7D 



refers to the paths from the exclusively defined class and the "knot entropy" problem arises. 



B. Abelian problems in statistics of entangled random walks and incompleteness of 

Gauss invariant 

As far back as 1967 S.F. Edwards had discovered the basis of the statistical theory of 
entanglements in physical systems. In he proposed the way of exact calculating the par- 
tition function of self-intersecting random walk topologically interacting with the infinitely 
long uncrossible string (in 3D case) or obstacle (in 2D-case). That problem had been con- 
sidered in mathematical literature even earlier — see the paper |T0| for instance — but S.F. 



Edwards was apparently the first to recognize the deep analogy between abelian topological 
problems in statistical mechanics of the Markov chains and quantum-mechanical problems 
(like Bohm-Aharonov) of the particles in the magnetic fields. The review of classical results 
is given in |T^, whereas some modern advantages are discussed in Wl . 



The 2D version of the Edwards' model is formulated as follows. Take a plane with an 
excluded origin, producing the topological constraint for the random walk of length L with 
the initial and final points Tq and respectively. Let trajectory make n turns around the 
origin (fig.@). The question is in calculating the distribution function 'P„(ro, r^, L). 
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FIG. 2. Random walk on the plane near the single obstacle. 



In the said model the topological state of the path C is fully characterized by number 
of turns of the path around the origin. The corresponding abelian topological invariant is 
known as Gauss linking number and when represented in the contour integral form, reads 



Inv{r(s)} = G{C} = f 

Jc 



ydx — xdy 



where 



AW 



c + y"^ 



ex 



/ A(r)dr = 2nn + {} 

Jc 



e= (0,0,1) 



(11.8) 



(11.9) 



and is the angle distance between ends of the random walk. 

Substituting Eg. ( [11. 8] ) into Eg. ( [11. 6] ) and using the Fourier transform of the 5-function, 
we arrive at 



P„(ro,ri,L) 



nLa 



exp 



La 




(11.10) 



which reproduces the well known old result (some very important generalizations one can 



find in 11 



Physically significant guantity obtained on the basis of Eg.( |11.10| ) is the entropic force 

d 



fnip) 



dp 



In P„(p,L) 



(11.11) 



which acts on the closed chain (fq = = p, ^ = 0) when the distance between the obstacle 
and a certain point of the trajectory changes. Apparently the topological constraint leads 
to the strong attraction of the path to the obstacle for any n ^ and to the weak repulsion 
for n = 0. 

Another exactly solvable 2D-problem closely related to the one under discussion deals 
with the calculation of the partition function of a random walk with given algebraic area. 
The problem concerns the determination of the distribution function Vs{ro,rL, L) for the 
random walk with the fixed ends and specific algebraic area 5*. 

As a possible solution of that problem, D.S. Khandekar and F.W. Wiegel |T3[ again repre- 
sented the distribution function in terms of the path integral Eg. (|II.6|) with the replacement 



(5[Inv{r(s)} - Inv] ^ 5[5{r(s)} - S] 



(11.12) 
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where the area is written in the Landau gauge: 

5'{r(s)} = - f ydx — xdy = - [ A{r}rds] A = ^ x r (11.13) 

(compare to Eqs.(jL8|)-(|L9D). 

The final expression for the distribution function reads ( |12[) 



where 



1 

Vsiro, TL, L) = — dg e^"^ P,(ro, rz., L) (11.14) 

ZTT J-oo 



(11.15) 



Pg(ro,rL,L) = 

A-K sm ^ 

X exp |^(xoyL - VoXl) - ^ ((xl - xof + {yi - y^f) cot 
and A = —iq. 

For closed trajectories Eqs. (p.l4|) -( [1I.15|) can be simplified essentially, giving 

ViiN) = 2La cosh^ (^^^ (11.16) 



Different aspects of this problem have been extensively studied in [11 



There is no principal difference between the problems of random walk statistics in the 
presence of a single topological obstacle or with a fixed algebraic area — both of them have 
the "abelian" nature. Nevertheless we would like to concentrate on the last problem because 
of its deep connection with the famous Harper-Hofstadter model dealing with spectral prop- 
erties of the 2D electron hopping on the discrete lattice in the constant magnetic field |]14 



Actually, rewrite Eq.( [11.4| ) with the substitution Eq.( [11.12D in form of recursion relation in 
the number of steps, N: 



(11.17) 



^^(rjv+i, N + 1) = J dvNQ (fat+i - r^) exp (-^^{^n x r/ 

For the discrete random walk on we use the identity 

/ dr^g {tn+i - )(...) ^ ^ w (ftv+i - r^) {■ ■ ■) (11.18) 

{riv} 

where w (r^v+i — tat) is the matrix of the local jumps on the square lattice; w is supposed 
to be symmetric: 

f I for (x, y) ^ {x,y± 1) and (x, y) ^ {x ± 1, y) 
w = I (11.19) 
[ otherwise 

Finally, we get in the Landau gauge: 
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-W{x,y,q,e) = e^''^''W{x,y - l,q) + e~^''^''W{x,y + l,q) + 

^ (11.20) 

e¥iyw{x - 1, y, q) + e-5^'?W(x + 1, y, q) 
where Vr(x,y,g, e) is the generating function defined via relation 

oo 

W{x,y,q,e)= e'^Vsivr,, N) 

N=0 

and q plays a role of the magnetic fiux through the contour bounded by the random walk 
on the lattice. 

There is one point which is still out of our complete understanding. On the one hand 
the continuous version of the described problem has very clear abelian background due to 
the use of commutative "invariants" like algebraic area Eq.( |ll.l3D . On the other hand it 
has been recently discovered ( |jl5|) that so-called Harper equation, i.e. Eq. (|II.20|) written 



in the gauge S{r} = J^ydx, exhibits the hidden quantum group symmetry related to the 
so-called C*-algebra ( |T6[) which is strongly nonabelian. Usually in statistical physics we 



expect that the continuous limit (when lattice spacing tends to zero with corresponding 
rescaling of parameters of the model) of any discrete problem does not change the observed 
physical picture, at least qualitatively. But for the considered model the spectral properties 
of the problem are extremely sensitive to the actual physical scale of the system and depend 
strongly on the lattice geometry. 

The generalization of the above stated problems concerns, for instance, the computation 
of the partition function for the random walk entangled with k > 1 obstacles on the plane 
located in the points {ri, . . . , Vk}- At first sight, approach based on usage of Gauss linking 
number as topological invariant, might allow us to solve such problem easily. Let us replace 
the vector potential A(r) in Eq.( [11.8|) by the following one 



A(r„...,r,)=ex^^-^ (11.21) 

j=i 1^ M1 

The topological invariant in this case will be the algebraic sum of turns around obstacles, 
which seems to be a natural generalization of the Gauss linking number to the case of 
many-obstacle entanglements. 

However, the following problem is bound to arise: for the system with two or more 
obstacles it is possible to imagine closed trajectories entangled with a few obstacles together 
but not entangled with every one. In the fig.^ the so-called "Pochhammer contour" is 
shown. Its topological state with respect to the obstacles cannot be described using any 
abelian version of the Gauss-like invariants. 




FIG. 3. Pochammer contour entangled with two obstacles together but not entangled with every 
one. 
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To clarify the point we can apply to the concept of the homotopy group |T^. Consider 
the topological space TZ = — {ri,r2} where {ri,r2} are the coordinates of the removed 
points (obstacles) and choose an arbitrary reference point Fq. Consider the ensemble of all 
directed trajectories starting and finishing in the point tq. Take the basis loops 7i(s) and 
72 (s) (0 < s < L) representing the right-clock turns around the points ri and r2 respectively. 
The same trajectories passed in the counter-clock direction are denoted by 7]~^(s) and 72^^ (s). 

The multiplication of the paths is their composition: for instance, 7172 = 71 o 72- The 
unit (trivial) path is the composition of an arbitrary loop with its inverse: 

e = 7.7-' = 7-S ^ = {1,2}. (11.22) 

The loops 7i(s) and 7j(s) are called equivalent if one can be transformed into another by 
means of monotonic change of variables s = s{s). The homotopic classes of directed trajec- 
tories form the group with respect to the paths multiplication; the unity is the homotopic 
class of the trivial paths. This group is known as the homotopy group 7ri(7?-, fq). 

Any closed path on 71 can be represented by the "word" consisting of set of letters 
{71, 72, 7]~^, 7^^}. Taking into account Eg. ( |11.22| ), we can reduce each word to the minimal 
irreducible representation. For example, the word W = 7i72'^7i7i7r^72'"^727r"^72'^ can be 
transformed to the irreducible form: W = 7i72'^72'^- K is easy to understand that the word 
W = e represents only the unentangled contours. The entanglement in fig.^ corresponds to 
the irreducible word W = 7r^727i72^^ = 1- The non-abelian character of the topological 
constraints is reflected in the fact that different entanglements do not commute: 7172 7^ 7271- 
At the same time, the total algebraic number of turns (Gauss linking number) for the 
path in fig.|^ is equal to zero, i.e. it belongs to the trivial class of homology. Speaking 
more formally, the mentioned example is the direct consequence of the well known fact in 
topology: the classes of homology of knots (of entanglements) do not coincide in general with 
the corresponding homotopic classes. The first ones for the group tti can be distinguished 
by the Gauss invariant, while the problem of characterizing the homotopy class of a knot 
(entanglement) by an analytically defined invariant is one of the main problems in topology. 

The principal difficulty connected with application of the Gauss invariant is due to its 
incompleteness. Hence, exploiting the abelian invariants for adequate classification of topo- 
logically different states in the systems with multiple topological constraints is very prob- 
lematic. 



C. Nonabelian Algebraic Knot Invariants 

The most obvious topological questions concerning the knotting probability during the 
random closure of the random walk cannot be answered using the Gauss invariant due to 
its weakness. 

The break through in that field was made in 1975-1976 when the algebraic polynomials 
were used for the topological state identification of closed random walks generated by the 
Monte-Carlo method |1[]. It has been recognized that the Alexander polynomials being much 
stronger invariants than the Gauss linking number, could serve as a convenient tool for the 
calculation of the thermodynamic properties of entangled random walks. That approach ac- 
tually appeared to be very fruitful and the main part of our modern knowledge on knots and 
links statistics is obtained with the help of these works and their subsequent modifications. 
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In the present Section we develop the analytic approach in statistical theory of knots 
considering the basic problem — the probability to find a randomly generated knot in a spe- 
cific topological state. We would like to reiterate that our investigation would be impossible 
without utilizing of algebraic knot invariants discovered recently. Below we reproduce briefiy 
the construction of Jones invariants following the Kauffman approach in the general outline. 



1. Disordered Potts model and generalized dichromatic polynomials 

The graph expansion for the Potts model with the disorder in the interaction constants 
can be defined by means of slight modification of the well known construction of the ordinary 
Potts model ||T8|,19|. Let us recall the necessary definitions. 



Take an arbitrary graph C with vertices. To each vertex of the given graph we 
attribute the "spin" variable cxj {i G [1,A^]) which can take q states labelled as 1,2, 
on the simplex. Suppose that the interaction between spins belonging to the connected 
neighboring graph vertices only contributes to the energy. Define the energy of the spin's 
interaction as follows 

J \ \ Jki OTfc = dz, (dfc, a;) - neighbors 

Eki = Jki o{ak, ai) = < (11.23) 

1 otherwise 

where Jki is the interaction constant which varies for different graph edges and the equality 
cTfc = CT; means that the neighboring spins take equal values on the simplex. 
The partition function of the Potts model now reads 

Zpotts = I] exp J ^ ^S{ak, cTi) \ (11.24) 

where T is the temperature. 

Expression Eq. ([11.2^ gives for g = 2 the well-known representation of the Ising model 
with the disordered interactions extensively studied in the theory of spin glasses |^ . (Later 
on we would like to fill in this old story by a new "topological" sense.) 

To proceed with the graph expansion of the Potts model |T^, rewrite the partition 
function ( |II.24| ) in the following way 

Zpotts = E n [1 + ^fc' ^i^k, ai)] (11.25) 

W} {kl} 

where 

Vki = exp (^^) - 1 (11.26) 

If the graph C has edges then the product Eq. ([II.25|) contains multipliers. Each 
multiplier in that product consists of two terms {1 and VkiS{ak,(7i)}. Hence the partition 
function Eq. ( P^I.25| ) is decomposed in the sum of 2^ terms. 

Each term in the sum is in one-to-one correspondence with some part of the graph C 
To make this correspondence clearer, it should be that an arbitrary term in the considered 
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sum represents the product of multipliers described above in ones from each graph edge. 
We accept the following convention: 

(a) If for some edge the multiplier is equal to 1, we remove the corresponding edge from 
the graph £; 

(b) If the multiplier is equal to vm 6{ak, cri) we keep the edge in its place. 

After repeating the same procedure with all graph edges, we find the unique representa- 
tion for all terms in the sum Eq. ( |1I.25| ) by collecting the components (either connected or 
not) of the graph C. 

Take the typical graph G consisting of m edges and C connected components where 
the separated graph vertex is considered as one component. The presence of 5-functions 
ensures the spin's equivalence within one graph component. As a result after summation of 
all independent spins and of all possible graph decompositions we get the new expression 
for the partition function of the Potts system Eq. (P!I.24|) 

m 

Zpotts = E ^TI ^ki (11.27) 

{G} {kl} 

where the product runs over all edges in the fixed graph G. 

It should be noted that the graph expansion Eq.(|L2^) where vm = v for all {k,l} 



coincides with the well known representation of the Potts system in terms of dichromatic 



polynomial (see, for instance, |T8|JT 



Another comment concerns the number of spin states, q. As it can be seen, in the 
derivation presented above we did not account for the fact that q has to take positive 
integer values only. From this point of view the representation Eq. (|II.271) has an advantage 
with respect to the standard representation Eq.([11.24D and can be considered an analytic 



continuation of the Potts system to the non-integer and even complex values of q. We 
show in the subsequent sections how the defined model is connected to the algebraic knot 
invariants. 



2. Reidemeister Moves and State Model for Construction of Algebraic Invariants 

Let K he a. knot (or link) embedded in the 3D-space. First of all we project the knot (link) 
onto the plane and obtain the 2D-knot diagram in the so-called general position (denoted by 
K as well). It means that only the pair crossings can be in the points of paths intersections. 
Then for each crossing we define the passages, i.e. parts of the trajectory on the projection 
going "below" and "above" in accordance with its natural positions in the 3D-space. 

For the knot plane projection with defined passages the following theorem is valid: (Rei- 



demeister [|22|): 

Two knots embedded in can be deformed continuously one into the other if and only if 
the diagram of one knot can be transformed into the diagram corresponding to another knot 
via the sequence of simple local moves of types I, II and III shown in fig.^j. 
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FIG. 4. Reidemeister moves of types I, II and III. 



The work provides us with the proof of this theorem. Two knots are called regular 
isotopic if they are isotopic with respect to two last Reidemeister moves (II and III); mean- 
while, if they are isotopic with respect to all moves, they are called ambient isotopic. As it 
can be seen from fig.^ the Reidemeister move of type I leads to the cusp creation on the 
projection. At the same time it should be noted that all real 3D-knots (links) are of ambient 
isotopy. 

Now, after the Reidemeister theorem has been formulated, it is possible to describe the 
construction of polynomial "bracket" invariant in the way proposed by L.H. Kauffman |7|,|23 



This invariant can be introduced as a certain partition function being the sum over the set 
of some formal ( "ghost" ) degrees of freedom. 

Let us consider the 2D-knot diagram with defined passages as a certain irregular lattice 
(graph). Crossings of path on the projection are the lattice vertices. Turn all these crossings 
to the standard positions where parts of the trajectories in each graph vertex are normal to 
each other and form the angles of ±7r/4 with the x-axis. It can be proven that the result 
does not depend on such standardization. 

There are two types of vertices in our lattice — a) and b) which we label by the variable 
bi = ±1 as it is shown below: 




h = +l 



and 



(b) 




The next step in the construction of algebraic invariant is introduction of two possible 
ways of vertex splittings. Namely, we attribute to each way of graph splitting the following 
statistical weights: A to the horizontal splitting and B to the vertical one for the vertex of 
type a); -B to the horizontal splitting and A to the vertical one for the vertex of type b). 
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The said can be schematically reproduced in the following picture: 

A ^ B 





^ A 

the constants A and B to be defined later. 

For the knot diagram with vertices there are 2^ different microstates, each of them 
representing the set of splittings of all vertices. The entire microstate, S", corresponds to 
the knot (link) disintegration to the system of disjoint and non-selfintersecting circles. The 
number of such circles for the given microstate S we denote as S. The following statement 
belongs to L. Kauffman ( [0)- 

Consider the partition function 

{K) = Y,(f-^A'B\ (11.28) 

{S} 

where means summation over all possible 2^ graph splittings, i and j = N — i being 

the numbers of vertices with weights A and B for the given realization of all splittings in the 
microstate S respectively. 

The polynomial in A, B and d represented by the partition function Eg. ( [11.2^ ) is the 
topological invariant of knots of regular isotopy if and only if the following relations among 
the weights A, B and d are fulfilled: 

AB = 1 

(11.29) 

ABd + A^ + B^ = 

The sketch of the proof is as follows. Denote with (. . .) the statistical weight of the knot 
or of its part. The (A')-value equals the product of all weights of knot parts. Using the 
definition of vertex splittings, it is easy to test the following identities valid for unoriented 
knot diagrams 



(11.30) 





completed by the "initial condition" 

'K[jO) = d{K); a: is not empty (11.31) 



where O denotes the separated trivial loop. 
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The skein relations Eq. (p.30|) correspond to the above defined weights of horizontal 
and vertical splittings while the relation Eg. ( [11.31] ) defines the statistical weights of the 
composition of an arbitrary knot and a single trivial ring. These diagrammatic rules are 
well defined only for fixed "boundary condition" of the knot (i.e., for the fixed part of the 
knot outside the brackets). Suppose that by convention the polynomial of the trivial ring is 
equal to the unity; 

O) = 1 (11.32) 

Now it can be shown that under the appropriate choice of the relations between A, B and 
d, the partition function Eq. ([11.28|) represents the algebraic invariant of the knot. The proof 



is based on direct testing of the invariance of (i^')-value with respect to the Reidemeister 
moves of types II and III. For instance, for the Reidemeister move of type II we have: 




52 



{ABd + A^ + B''' 




AB 



(11.33) 



Therefore, the invariance with respect to the Reidemeister move of type II can be ob- 
tained immediately if we set the statistical weights in the last line of Eq. ( |11.33D as it is 
written in Eq. ( [11.29]) . Actually, the topological equivalence of two knot diagrams is restored 
with respect to the Reidemeister move of type II only if the right- and left-hand sides of 
Eq.( [11.33[) are identical. It can also be tested that the condition of obligatory invariance 
with respect to the Reidemeister move of type III does not violate the relations Eq.( [il.29[ ). 

The relations Eq. ([11.29|) can be converted into the form 



B = A- 



d 



-A' - A- 



(11.34) 



which means that the Kauffman invariant Eq.( [11.28[ ) is the Laurent polynomial in A- value 
only. 

Finally, Kauffman showed that for oriented knots (links) the invariant of ambient isotopy 
(i.e., the invariant with respect to all Reidemeister moves) is defined via relation: 



m] = i-A) 



3Tw{K) 



(11.35) 



here Tw{K) is the twisting of the knot (link), i.e. the sum of signs of all crossings defined 
by the convention: 
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(not to be confused with the definition of the variable bi introduced above). Eg. ( [11.35]) 
follows from the following chain of equalities 



/ 





The state model and bracket polynomials introduced by L.H. Kauffman seem to be very 
special. They explore only the peculiar geometrical rules such as summation over the formal 
"ghost" degrees of freedom — all possible knot (link) splittings with simple defined weights. 
But one of the main advantages of the described construction is connected with the fact 
that Kauffman polynomials in A-value coincide with Jones knot invariants in t-value (where 
t = A'/'). 

Jones polynomial knot invariants were discovered first by V.F.R. Jones during his in- 
vestigation of topological properties of braids (see Section 3 for details). Jones' proposition 
concerns the establishment of the deep connection between the braid group relations and 
the Yang-Baxter equations ensuring the necessary condition of transfer matrix commutativ- 
ity 0. The Yang-Baxted equations play an exceptionally important role in the statistical 
physics of integrable systems (such as ice, Potts, 0{n), 8- vertex, quantum Heisenberg models 

The attempt to apply Kauffman invariants of regular isotopy for investigation of statis- 
tical properties of random walks with topological constraints in a thin slit has been made 



recently [24|. Below we extend the ideas of the work |24| considering the topological state 



of the knot special kind of a quenched disorder. 



D. Lattice knot diagrams as disordered Potts model 

Let us specify the model under consideration. Take a square lattice turned to the 
angle 7r/4 with respect to the x-axis and project a knot embedded in onto Ai supposing 
that each crossing point of the knot diagram coincides with one lattice vertex without fall 
(there are no empty lattice vertices) — see fig.|^. Define the passages in all N vertices and 
choose such boundary conditions which ensure the lattice to form a single closed path; that 
is possible when \/N (i.e. A^) is an odd number. The frozen pattern of all passages {bi} on 
the lattice together with the boundary conditions fully determine the topology of some 3D 
knot. 
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FIG. 5. Lattice knot with topological disorder realized in a quenched random pattern of passages. 

Of course, the model under consideration is rather rough because we neglect the "space" 
degrees of freedom due to trajectory fluctuations and keep the pure topological specificity 
of the system. Later on in Chapter 4 we discuss the applicability of such model for real 
physical systems and produce arguments in support of its validity. 

The basic question of our interest is as follows: what is the probability VN{f[K]} to 
find a knot diagram on our lattice in a topological state characterized by some specific 
Kauffman invariant f[K] among all 2^ microrealizations of the disorder in the lattice 
vertices. That probability distribution reads (compare to Eq.( [ll.l|) ) 

1 



(11.36) 



{b^} 



where f[K{bi, . . . , ^at}] is the representation of the Kauffman invariant as a function of all 
passages {bi} on the lattice A4. These passages can be regarded as a sort of quenched 
"external field" (see below). 

Our main idea of dealing with Eq.( [lI.36D consists in two steps: 

(a) At first we convert the Kauffman topological invariant into the known and well- 
investigated Potts spin system with the disorder in interaction constants; 

(b) Then we apply the methods of the physics of disordered systems to the calculation of 
thermodynamic properties of the Potts model. It enables us to extract finally the estimation 
for the requested distribution function. 

Strictly speaking, we could have disregarded point (a), because it does not lead directly 
to the answer to our main problem. Nevertheless we follow the mentioned sequence of 
steps in pursuit of two goals: 1) we would like to prove that the topologically-probabilistic 
problem can be solved within the framework of standard thermodynamic formalism; 2) we 
would like to employ the knowledges accumulated already in physics of disordered Potts 
systems to avoid some unnecessary complications. Let us emphasize that the mean-field 
approximation and formal replacement of the model with short-range interactions by the 
model with infinite long-range ones serves to be a common computational tool in the theory 
of disordered systems and spin glasses. 
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1. Algebraic invariants of regular isotopy 



The general outline of topological invariants construction deals with seeking for the 
functional, f[K{bi,...bN}], which is independent on the knot shape i.e. is invariant with 
respect to all Reidemeister moves. 

Recall that the Potts representation of the Kauffman polynomial invariant Eq. ( [1I.28| ) of 
regular isotopy for some given pattern of "topological disorder", {bi}, deals with simulta- 
neous splittings in all lattice vertices representing the polygon decomposition of the lattice 
Ai. Such lattice disintegration looks like a densely packed system of disjoint and non- 
selfintersecting circles. The collection of all polygons (circles) can be interpreted as a system 
of the so-called Eulerian circuits completely filling the square lattice. Eulerian circuits are 
in one-to-one correspondence with the graph expansion of some disordered Potts system 
introduced in Section 2.3.1 (see details below and in [p7[]). 

Rewrite the Kauffman invariant of regular isotopy, {K) , in form of disordered Potts model 
defined in the previous section. Introduce the two-state "ghost" spin variables, = ±1 in 
each lattice vertex independent on the crossing in the same vertex 




Irrespective of the orientation of the knot diagram shown in fig.|^ (i.e. restricting with 
the case of regular isotopic knots), we have 

{K{h}) = E + ^"')^"'exp [\nAj2hs\ (11.37) 

{S} V i=l I 

Written in such form the partition function represents the weihgted sum of all 

possible Eulerian circuits^ on the lattice M.. Let us show explicitly that the microstates 
of the Kauffman system are in one-to-one correspondence with the microstates of some 
disordered Potts model on a lattice. Apparently for the first time the similar statement was 
expressed in the paper 0]. To be careful, we would like to use the following definitions: 

(i) Let us introduce the lattice L dual to the lattice M., or more precisely, one of two 
possible (odd and even) diagonal dual lattices, shown in fig.||. It can be easily noticed that 
the edges of the lattice L are in one-to-one correspondence with the vertices of the lattice 
M.. Thus, the disorder on the dual lattice L is determined on the edges. In turn, the edges 
of the lattice L can be divided into the subgroups of vertical and horisontal bonds. Each 
fcZ-bond of the lattice L carries the "disorder variable" b^i being a function of the variable bi 
located in the corresponding i- vertex of the lattice M.. The simplest and most sutable choice 
of the function buibi) is as in Eq. ( [11.4^ ) (or vice versa for another choice of dual lattice); i 
is the vertex of the lattice M. belonging to the fc/-bond of the dual lattice L. 



Eulerian circuit is a trajectory on the graph which visits once and only once all graph edges. 
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polygon decomposition 
of lattice J4 



backbone graph 
on dual lattice £ 




FIG. 6. Disintegration of the knot diagram on the 7W -lattice into ensemble of nonselfintersecting 
loops (Eulerian circuits) and graph representation of the Potts model on the dual £-lattice. 

(ii) For the given configuration of splittings on M. and chosen dual lattice £ let us accept 
the following convention: we mark the edge of the £-lattice by the solid line if this edge 
is not intersected by some polygon on the A^-latice and we leave the corresponding edge 
unmarked if it is intersected by any polygon — as it is shown in the fig.^. Similarly, the 
sum J^Sibi in Eq.( |11.37| ) can be rewritten in terms of marked and unmarked bonds on the 
^-lattice 



Sibi + ^i^i 

mark nonmark 

horiz vertic horiz vertic 

mark mark nonmark nonmark 

horiz vertic horiz vertic 

- Y ^kl - Y + Y ^kl + Y ^kl 

mark mark nonmark nonmark 



(11.38) 



Y bki-Y = Y 



nonmark 



mark 



bki — 2 Y, bki 

all edges mark 



where we used the relation J2 b^i + J2 bki = J2 b^i. 

nonmark mark all edges 

(iii) Let rris be the number of marked edges and Cs be the number of connected compo- 
nents of marked graph. Then the Euler relation reads: 



S = 2Cs + ms~N + x 



(11.39) 



The Eq.( [11.39 ) can be proved directly. The X" value depends on the genus of the surface, 
which can be covered by the given lattice, (i.e. x depends on the boundary conditions). In 
the thermodynamic limit ^ 1 the x-dependence should disappear (at least for the flat 
surfaces), so the standard equality S = 2Cs + mg — N will be assumed below. 



20 



By means of definitions (i)-(iii), we can easily convert Eq. (p!I.37|) into the form: 



TV 



all edges 

nis 



(11.40) 



{G} 



mark 



where we used Eg. ( |11.38D and the fact that + 1 is even. Comparing Eq.( [11.40| ) with 
Eg. ( [II. 25]) we immediately conclude that 



Ills 

J2{A^ + A-'f""^ n [A-^'^'i-A'-A-^] 



= J2Ili'^ + ^kiS{cxk,cxi)) (11.41) 

{a} {kl} 



what coincides with the partition function of the Potts system written in the form of dichro- 
matic polynomial. Therefore, we have 



Vki=^A-^'^'(-A^-A-^) 



.1 _ y4-4ftfc; 



q={A^ + A 



-2\2 



(11.42) 



Since the "disorder" variables bki take the discrete values ±1 only, we get the following 
expression for the interaction constant Jki (see Eg. ( [11.26]) ) 



Jkl 

T 



In 



1-{A^ + A-^)A 



\n[-A-^''^'] 



(11.43) 



Combining Egs.( 11.40 )-( 11.43|) we obtain the following statement. 

(a) Take N -vertex knot diagram on the lattice M. with given boundary conditions and 
fixed set of passages {hi}, (b) Take the dual lattice C in one-to-one correspondence with M. 
where one vertex of M. belongs to one edge of C 

The Kaufjman topological invariant {K{A)) of regular isotopy for knot diagrams on Ai 
admits representation in form of 2D Potts system on the dual lattice C: 



where: 



{K{A)) = H{A,{hM}) ZpoUs[q{A),{MhkuA)}] 



H{A, {bki}) = {A^ + A-^Y^^^''' exp | In A ^ b^ 

{kl} 



(11.44) 



(11.45) 



is the trivial multiplier [H does not depend on Potts spins); 



Zpotts[q{A), {Jki{bki, A)}] = ^exp IJ2 ^^^%r— ^'^('^fc^c^O 



.{kl} 



T 



(11.46) 



is the Potts partition function with interaction constants, Jm, and number of spin states, q, 
defined as follows 



21 



^ = ln[-/l-^^-]; q={A^ + A'^)' (11.47) 

and the variables bki play a role of disorder on edges of the lattice C dual to the lattice A4 . 
The connection between bki and bi is defined by convention 

{—bi if (kl)- edge is vertical 
(11.48) 
bi if (kl)- edge is horizontal 

Eq.( |11.4l|) has the sense of partition function of the 2D disordered Potts system with 
the random nearest-neighbor interactions whose distribution remains arbitrary. The set of 
passages {bki} uniquely determines the actual topological state of the woven carpet for the 
definite boundary conditions. Therefore the topological problem of the knot invariant deter- 
mination is reduced to usual statistical problem of calculation of the partition function of the 
Potts model with the disorder in the interaction constants. Of course, this correspondence 
is still rather formal because the polynomial variable A is absolutely arbitrary and can take 
even complex values, but for some regions of A that thermodynamic analogy makes sense 
and could be useful as we shall see below. 

The specific feature of the Potts partition function which gives the representation of the 
Kauffman algebraic invariant is connected with the existence of the relation between the 
temperature T and the number of spin states q (see Eq.( [11.42| )) according to which T and q 
cannot be considered anymore as independent variables. 



2. Algebraic invariants of ambient isotopy 

The invariance of the algebraic topological invariant, f[K], with respect to all Reide- 
meister moves (see Eq.(|L35)) for our system shown in the fig.^ is related to the oriented 



Eulerian circuits called Hamiltonian walks 0. 

Let us suppose that the orientation of the knot diagram shown in fig.|^ is chosen according 
to the natural orientation of the path representing a knot K mP? . For the defined boundary 
conditions we get the so-called Manhattan lattice consisting of woven threads with alternating 
directions. 

It follows from the definition of twisting Tw{K) (see the Section 1.3.2) that Tw{K) 
changes the sign if the direction of one arrow in the vertex is changed to the inverse. Revers- 
ing the direction of any arrows in the given vertex even times we return the sign of twisting 
to the initial value. 

We define groups of "even" and "odd" vertices on the lattice Ad as follows. The vertex 
i is called even (odd) if it belongs to the horizontal (vertical) bond (kl) of the dual lattice 
C. Now it is easy to prove that the twisting of the knot on the Manhattan lattice Ai can 
be written in terms of above defined variables bki- Finally the expression for the algebraic 
invariant of ambient isotopy f[K] on the lattice C reads 



Hamiltonian walk is a closed path which visits once and only once all vertices of the given 
oriented graph. 
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f[K]=exp\3\n[-A]Y,hi] {K{{bM},A)) (11.49) 

m 



where {K{{bki},A)) is defined by Eq.([mD 



E. Notion about annealed and quenched realizations of topological disorder 

Fixed topological structure of a trajectory of given length fiuctuating in space is a typical 
example of a quenched disorder. Actually, the knot structure is formed during the random 
closure of the path and cannot be changed without the path rupture. Because of the topo- 
logical constraints the entire phase space of ensemble of randomly generated closed loops is 
divided into the separated domains resembling the multi-valley structure of the spin glass 
phase space. Every domain corresponds to the sub-space of the path configurations with 
the fixed value of the topological invariant. 

The methods of theoretical description of the systems with quenched disorder in inter- 
action constants are rather well developed, especially in regard to the investigation of spin 



glass models |20]. Central for these methods is the concept of self- averaging which can be 
explained as follows. Take some additive function F (the free energy, for instance) of some 
disordered spin system. The function F is the self-averaging quantity if the observed value. 
Fobs, of any macroscopic sample of the system coincides with the value Fav averaged over 
the ensemble of disorder realizations: 



obs 



The central technical problem is in calculation of the free energy F = —T In Z averaged over 
the randomly distributed quenched pattern in the interaction constants. In this Section we 
show that this famous thermodynamic problem of the spin glass physics is closely related to 
the knot entropy computation. 

Another problem arises when averaging the partition function Z (but not the free energy) 
over the disorder. Such problem is much simpler from computational point of view and 
corresponds to the case of annealed disorder. Physically such model corresponds to the 
situation when the topology of the closed loop can be changed. It means that the topological 
invariant, i.e. the Potts partition function, has to be averaged over all possible realizations 
of the pattern disorder in the ensemble of open (i.e. unclosed) loops on the lattice. It has 
been shown in that the calculation of the mean values of topological invariants allows 
to extract rather rough but nontrivial information about the knot statistics. 



1. Entropy of knots. Replica methods 

Our main goal is the computation of the probability distribution VN{f[K]} (see 
Eq.( [lI.36D ). Although we are unable to evaluate this function exactly, the representation 
of VN{f[K]} in terms of disordered Potts system enable us to give an upper estimation 
for the fraction of randomly generated paths belonging to some definite topological class (in 
particular, to the trivial one). We use the following chain of inequalities restricting ourselves 
with the case of regular isotopic knots for simplicity ( |]2^): 
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Probability Vn of 
knot formation in 
a given topological 

state 



< 



Probability Vn{K{A)} of 
knot formation with specific 
topological invariant {K{A)) 
for all A 



< 



Probability Vn{K{A*)} of 
knot formation for specific 
value of A* minimizing the 
free energy of associated 
Potts system 



(11.50) 

The first inequality is due to the fact that Kauffman invariant of regular isotopic knots 
is not a complete topological invariant, whereas the last probability in the chain can be 
written as follows 



Vn{K{A*)} = ^ Q{bM}A[{K{hi,A*}) - {K{A*)) 

{bkl} 



(11.51) 



where J2 means summation over all possible configurations of the "crossing field" {bki}, 
A-function cuts out all states of the field {bki} with specific value of Kauffman invariant 
{K{bki, A*}) = {K{A*)) and 0{6h} is the probability of realization of given crossings con- 
figuration. 

In principle the distribution 6(6,} depends on statistics of the path in underlying 3D 
space and is determined physically by the process of the knot formation. Here we restrict 
ourselves to the following simplest suppositions: 

(i) We regard crossings {bi} in different vertices of A^-lattice as completely uncorrelated 
variables (or, in other words, we assume that the variables {bki} defined on the edges of the 
^-lattice are statistically independent): 



N 



(11.52) 



(ii) We suppose variable bi (or bki) to take values ±1 with equal probabilities, i.e.: 



Pih) = \ Sih + l 5{h + 1) 
The probability of trivial knot formation can be estimated now as follows 
rP{A*) < J2'S>{bki}A\\n{K{bki,A*}) 



(11.53) 



{bkl} 



(11.54) 



EPibM)dbki {K^nhi,A*} 

kl 



where {K{A*)) = 1 for trivial knots. 

Thus our problem is reduced to the calculation of non-integer complex moments of the 
partition function, i.e., values of the type {K^y{bki, A*}). An analogous problem of evaluation 
of non-integer moments is well known in the spin-glass theory. Indeed, the averaging of the 
free energy of the system, F, over quenched random field is widely performed via so-called 

nlnZ 



replica-trick p8[. The idea of the method is as follows. Consider the identity = e 
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and expand the right-hand side up to the first order in n. We get = 1 + n\n Z + 0{n'^). 
Now we can write 

Z" — 1 

F = —\nZ = — hm 

n^O n 

We proceed with the calculation of the complex moments of the partition function 
{K{bki})- In other words we would like to find the averaged value (K^) for integer val- 
ues of n. Then we put n = iy and compute the remaining integral in Eg. ( |ll-5^ ) over value. 
Of course, this procedure needs to be verified and it would be of most desire to compare our 
predictions with the data obtained in numerical simulations. However let us stress that our 
approach is no more curious than replica one, it would be extremely desirable to test the 
results obtained by means of computer simulations. 

The outline of our calculations is as follows. We begin by rewriting the averaged Kauff- 
man invariant using the standard representation of the replicated Potts partition function 
and extract the corresponding free energy F{A) in the frameworks of the infinite-range 
mean-field theory in two dimensions. Minimizing F{A) with respect to A we find the equi- 
librium value A*. Then we compute the desired probability of trivial knot formation {A*) 
evaluating the remaining Gaussian integrals. 

Averaging the nth power of Kauffman invariant over independent values of the "crossing 
field" bki = ±1 we get 



l[P{bM)dbMK^''{bM} 

kl 



[2 cosh(2/?)] 



-2?i(Af+l) 



(11.55) 



E n E ^(^fc ' + In cosh 

{a} kl [ kl 



a=l 



where f3 = In A. Let us break for a moment the connection between the number of spin states, 
q, and interaction constant and suppose -C 1. Later on we shall verify the selfconsistency 
of this approximation. Now the exponent in the last expression can be expanded as a power 
series in f3. Keeping the terms of order only, we rewrite Eg. ( [[I.55D in the standard form 
of n-replica Potts partition function 



{K^{A)) = [2cosh(2/3)]-'"('^+')exp 



N ( Ip'n' 



J' 



2 N n 



E exp — 



^ka^kb^la^lb 



/ j2 \ N n 

h^(g-2) + JoK:E<^f. 

\ ^ J kl a=l 



(11.56) 



{o-i...(7„} kl a^f3 

where spin indexes a, b change in the interval [0, g — 1], ^ 1 and 

Jq = in — AjS'^n 
g = 4 + 16/^2 > 4 



(11.57) 
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According to the results of Cwilich and Kirkpatrick [|29| and later works (see, for in- 
stance, [^), the spin-glass ordering takes place and the usual ferromagnetic phase makes 
no essential contribution to the free energy under the condition 



Jo g-4 
J 2 



(11.58) 



Substituting Eq. (p37[) into Eq.(|L58D it can be seen that 3?(l.h.s.) < 3fJ(r.h.s) in Eq.( p38D 
for all j3. Thus, we expect that the spin-glass ordering (in the infinite-range model) corre- 
sponds to the solutions 



Q 



al3 
ab 



q6iala)-l) = 
q6K,a)-l){q5i4,b)-l)^0 



where m° and are the ferromagnetic and spin-glass order parameters respectively. If 
it is so, we can keep the term in the exponent (Eq.(|L5^)) corresponding to inter-replica 
interactions only. 

We follow now the standard scheme of analysis of Potts spin glasses partition function 



exhaustively described in p9| , |3ll , |30| ; main steps of this analysis are shortly represented below. 
Performing the Hubbard- Stratonovich transformation to the scalar fields and implying 



the homogeneous isotropic solution of the form Q^'^ 
JK^ (Eq. ( [11560 ) as follows ( g 



Qt^^ab, we can write down the value 



(iT") = exp<^ 



TT 

In —n{n 



l){q - 1)2 - In 2 cosh 



J 



~32 



ElUdQfexpLjHiQfjd' 

{a} i 



X 



(11.59) 



where 



H{<r"} = (q-l] 



q 



1 ^2 



(11.60) 



8 



g2 — 6g + 6 
48 



EW 

a/3 



In ||3T1 , [29[1 it was shown, that the mean-field replica symmetric solution of the mean-field 
Potts spin glass is unstable for g > 2 and the right ansatz of Eqs. (|II.59|) - (|II.60|) corresponds 
to the first level of Parisi replica breaking scheme for spin glasses. Hence, we have 



Q 



a/3 



Q if a and (3 belong to the same group of m replicas 



otherwise 



[n.61) 
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Analysis shows that for g > 4 (our case) the transition to the glassy state corresponds 
to m = 1 which implies the accessory condition Fpm = Fgg, where Fpm and Fgg are the free 
energies of paramagnetic and spin-glass phases respectively. The transition occurs at the 
point 

4)2 



1 - 



(11.62) 



J2 3(g2-18g + 42) 

Substituting Eg. ( |11.57l ) into Eq.( p.l.62]) we find the self-consistent value of reverse temperature 
of a spin-glass transition, Ptr- 



Ptr ~ 0.35 



(11.63) 



This numerical value is consistent with the condition <^ 1 implied above in the course of 
expansion of Eq. ( |11.56| ) . 

According to the results of the work the n-replica free energy near the transition 
point has the following form 

2 



64 



Nn{q-iyQ 



tr 



with the following expression of the spin-glass order parameter 

2(4 -g) 



Q 



tr 



> 



(11.64) 



(11.65) 



- 18q + 42 

From Eq. (|II.64|) we conclude that the free energy F reaches its minimum as a function 
of A = exp (P) just at the point A* = exp{Ptr). Using Eqs. (|11.64|) and (|II.65| ) we rewrite 
the expression for the averaged n-replica Kauffman invariant {K"^) in the vicinity of Ptr as 
follows (compare to ||29(|): 



(X2n^ ~ exp<^ Nn 



(3 + 16P' 



In 



TT 



I6/32 2 



Nn 



(3 + 16p' 



In 



71 , p^ 

i6^+^^' + y 



(3 + 16/32)^ (r^- a;') 'a 



(11.66) 



iW+42. 
iy, we can easily 



(4 + 16^2^^-18(4 + 16/3, 

Substituting Eq. ( [11.66]) into Eq. (|II.54|) and bearing in mind, that n 
evaluate the remaining Gaussian integral over value and obtain the result for V^^ {A) . As 
it has been mentioned above, to get the simplest estimation for probability of trivial knot 
formation, we use the last inequality in the chain of equations ( 11.50 ) corresponding to the 
choice A = A* = exp(/3jr): 



(11.67) 



This dependence it is not surprising from the point of view of statistical mechanics 
because the value 77 = V^\A*) is proportional to the free energy of the Potts system. But 
from the topological point of view the value rj has the sense of typical "complexity' of the 
knot (see also Section 3). The fact that t] grows linearly with means that the maximum 
of the distribution function Pif], N) is in the region of very "complex" knots, i.e. knots 
far from trivial. This circumstance directly follows from the non-commutative nature of 
topological interactions. 
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III. RANDOM WALKS ON LOCALLY NON-COMMUTATIVE GROUPS 



Recent years have been marked by the emergence of more and more problems related to 
the consideration of physical processes on non-commutative groups. In trying to classify such 
problems, we distinguish between the following categories in which the non-commutative 
origin of phenomena appear with perfect clarity: 

1. Problems connected with the spectral properties of the Harper-Hofstadter equation 
dealing with the electron dynamics on the lattice in a constant magnetic field. We mean 

primarily the consideration of groups of magnetic translations and properties of quantum 
planes [|2|,|T§. 

2. Problems of classical and quantum chaos on hyperbolic manifolds: spectral properties 
of dynamical systems and derivation of trace formulae as well as construction of 
probability measures for random walks on modular groups [Q. 

3. Problems giving rise to application of quantum group theory in physics: deformations 
of classical abelian objects such as harmonic oscillators 0] and standard random walks 0. 

4. Problems of knot theory and statistical topology: construction of nonabelian topo- 
logical invariants [p|,p3| , consideration of probabilistic behavior of the words on the simplest 
non-commutative groups related to topology (such as braid groups) ||3^, statistical proper- 
ties of "anyonic" systems [BS . 



5. Classical problems of random matrix and random operator theory and localization 
phenomena: determination of Lyapunov exponents for products of random non-commutative 



matrices p9|- ^] , study of the spectral properties and calculation of the density of states of 
large random matrices 

Certainly, such a division of problems into these categories is very speculative and reflects 
to a marked degree the authors' personal point of view. However, we believe that the 
enumerated items reflect, at least partially, the currently growing interest in theoretical 
physics of the ideas of non-commutative analysis. Let us stress that we do not touch upon 
the pure mathematical aspects of non- commutative analysis in this paper and the problems 
discussed in the present work mainly concern the points 4 and 5 of the list above. 

In the present Section we continue analyzing the statistical problems in knot theory, but 
our attention is paid to some more delicate matters related to investigation of correlations 
in knotted random paths caused by the topological constraints. The methods elaborated 
in Section 2 allow us to discuss these questions but we find it more reasonable to take a 
look at the problems of knot entropy estimation in terms of conventional random matrix 
theory. We believe that many non-trivial properties of the knot entropy problem can be 
clearly explained in context of the limit behavior of random walks over the elements of some 
non-commutative (hyperbolic) groups PB |. 

Another reason which forces us to consider the limit distributions (and conditional limit 
distributions) of Markov chains on locally non- commutative discrete groups is due to the 
fact that this class of problems could be regarded as the first step in a consistent harmonic 
analysis on the multiconnected manifolds (like Teichmiiller space); see also the Section 4. 
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A. Brownian bridges on simplest non-commutative groups and knot statistics 



As it follows from the said above the problems dealing with the investigation of the 
limit distributions of random walks on non-commutative groups is not a new subject in the 
probability theory and statistical physics. 

However in the context of "topologically-probabilistic" consideration the problems deal- 
ing with distributions of non- commutative random walks are practically out of discussion, 
except for very few special cases PT| , ^ , ^ . Particularly, in these works it has been shown 
that statistics of random walks with the fixed topological state with respect to the regular 
array of obstacles on the plane can be obtained from the limit distribution of the so-called 
"Brownian bridges" (see the definition below) on the universal covering — the graph with the 
topology of Cayley tree. The analytic construction of nonabelian topological invariant for 
the trajectories on the double punctured plane and statistics of simplest nontrivial random 
braid was shortly discussed in [El . 



Below we calculate the conditional limit distributions of the Brownian bridges on the 
braid group i?3 and derive the limit distribution of powers of Alexander polynomial of knots 
generated by random i?3-braids. We also discuss the limit distribution of random walks on 
locally free groups and express some conjectures about statistics of random walks on the 
group Bn. More extended discussion of the results concerning the statistics of Markov chains 



on the braid and locally free groups one can find in |52-54 



1. Basic definitions and statistical model 

The braid group i?„ of n strings has n — 1 generators {ai, (72, ... , with the following 

relations: 

cTiO-j+io-i = cTi+iO-jO-i+i {1 < i < n-1) 

aiCj = ajCi j\ > 2) (HI.l) 

cr-o-ri = o-rVi = e 

Any arbitrary word written in terms of "letters" — generators from the set {cti,..., 
o'n-ij crfS • • • 5 (^n-i} — gives a particular braid. The geometrical interpretation of braid gen- 
erators is shown below: 




1 2 ••• i i + 1 ■ ■ ■ n - 1 n 
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The length of the braid is the total number of the used letters, while the minimal 
irreducible length hereafter referred to as the "primitive word" is the shortest noncon- 
tractible length of a particular braid which remains after applying all possible group relations 
Eq.( [lII.l] ). Diagrammatically the braid can be represented as a set of crossed strings going 
from the top to the bottom appeared after subsequent gluing the braid generators. 

The closed braid appears after gluing the "upper" and the "lower" free ends of the braid 
on the cylinder. 

Any braid corresponds to some knot or link. So, it is feasible principal possibility to use 
the braid group representation for the construction of topological invariants of knots and 
links. However the correspondence between braids and knots is not mutually single valued 
and each knot or link can be represented by infinite series of different braids. This fact 
should be taken into account in course of knot invariant construction. 

Take a knot diagram K in general position on the plane. Let f[K] be the topological 
invariant of the knot K. One of the ways to construct the knot invariant using the braid 
group representation is as follows. 

1. Represent the knot by some braid b G -B„. Take the function / 

Demand / to take the same value for all braids b representing the given knot K. That 
condition is established in the well-known Markov-Birman theorem (see, for instance, [0): 
The function fxib} defined on the braid b E Bn is the topological invariant of a knot or 
link if and only if it satisfies the following "Markov condition" : 



fK{b'b"} = fK{b"b'} 

fxib' an] = fKWn b'} = fiAb'} b', b" e 5„ 
where b' and b" are two subsequent sub-words in the braid — see fig.^. 



(III.2) 



f 



= f 



f 



= f 



30 



FIG. 7. Geometric representation of Eqs.( |ll.2 ). 



2. Now the invariant fxib} can be constructed using the hnear functional if{b} defined 
on the braid group and called Markov trace. It has the following properties 

^{b'b"} = ip{b"b'} 

ifib' a^} = Tifib'} (III.3) 
V{b' a-'} = f^{b'} 

where 

T = ^{a,}, f = ^{cr7^}; ^e[l,n-l] (III.4) 

The invariant fxib} of the knot K is connected with the linear functional (p{b} defined 
on the braid b as follows 

Mb} = {rfr^^-'y^ (^-j ^{b} (III.5) 

where #(+) and #(— ) are the numbers of "positive" and "negative" crossings in the given 
braid correspondingly. 

The Alexander algebraic polynomials are the first well-known invariants of such type. 
In the beginning of 1980s Jones discovered the new knot invariants. He used the braid 
representation "passed through" the Hecke algebra relations, where the Hecke algebra, Hn{t), 
for Bn satisfies both braid group relations Eq. ([Ill.lj ) and an additional "reduction" relation 



(see the works |55| , |56| ) 

a^ = {l-t)ai + t (III.6) 

Now the trace (f{b} = (p{t){b} can be said to take the value in the ring of polynomials of one 
complex variable t. Consider the functional ip{t) over the braid {b' aib"}. Eq.( |111.6D allows 
us to get the recursion (skein) relations for ip{t) and for the invariant /x(^) (see for details 



ip{t){b'aib"} = (1 - t)if{t){b'b"} + t<f{t){b'(T-^b"} (111.7) 

and 

~ 1 ^2 

mt) - 1 iz) fKit) = a-t) (z) fUt) (iii.s) 



where = f{b'cxib"}; = f{b'cx-H"}; = f{b'b"} and the fraction - depends on 

T 

the used representation. 

3. The tensor representations of the braid generators can be written as follows 

a,iu) = Jim J2 R^inW^^^ ® ■ ■ ■ ^^'"'^ ® Kk ® ® ^^'^'^ ® ■ ■ ■ ^^"^ (III.9) 

klmn 
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where 7^*^ is the identity matrix acting in the position z; Enk is a matrix with {E, 



and Rf^ is the matrix satisfying the Yang-Baxter equation 



ka \ 



nkj'pq ^np^kq 



(III.IO) 



abc 



abc 



In that scheme both known polynomial invariants (Jones and Alexander) ought to be 
considered. In particular, it has been discovered in []57| , |58| that the solutions of Eq. (|III.10|) 
associated with the groups SUq{2) and GL{1, 1) are linked to Jones and Alexander invariants 
correspondingly. To be more specific: 

(a) — = for Jones invariants, fxif) = ^if)- The corresponding skein relations are 



r 



t-^v+{t) - tv-{t) = _ ^1/2)^0^^^ 



(iii.ii) 



and 



r 



are 



(b) — = t for Alexander invariants, fxit) = V(t). The corresponding skein relations 

r 



V+(t) - V"(t) = {t 



-1/2 



ti/2)V°(t) 



(III.12) 



To complete this brief review of construction of polynomial invariants from the repre- 
sentation of the braid groups it should be mentioned that the Alexander invariants allow 
also another useful description [^. Write the generators of the braid group in the so-called 
Magnus representation 



/I 

■•. 



\ 



A 



V 



■•. 

ly 



jth row; A 



'\ o\ 

t -t 1 



;iii.i3) 



Now the Alexander polynomial of the knot represented by the closed braid W = H^Li c"a 



of the length N one can write as follows 

{l+t + f + ... + e-^) V(t){A} = det 



N 

n ^c., - e 



(III.14) 



where index j runs "along the braid", i.e. labels the number of used generators, while the 
index a = {1, . . . , n — 1, n, . . . , 2n — 2} marks the set of braid generators (letters) ordered 
as follows {(Ti, . . . , cT„_i, (Tf ^, . . . , <Jn-i}- In our further investigations we repeatedly address 
to that representation. 



^Let us stress that the standard skein relations for Alexander polynomials one can obtain from 
Eg . (1111.121) replacing t^/"^ by -t'^/'^. 
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We are interested in the limit behavior of the knot or hnk invariants when the length 
of the corresponding braid tends to infinity, i.e. when the braid "grows". In this case we 
can rigorously define some topological characteristics, simpler than the algebraic invariant, 
which we call the knot complexity. 

Call the knot complexity, rj, the power of some algebraic invariant, fxit) (Alexander, 
Jones, HOMFLY) (see also f^) 

|t|-»oo ln|i;| 

Remark. By definition, the "knot complexity" takes one and the same value for rather 
broad class of topologically different knots corresponding to algebraic invariants of one and 
the same power, being from this point of view weaker topological characteristics than com- 
plete algebraic polynomial. Let us summarize the advantages of knot complexity. 

(i) One and the same value of rj characterizes a narrow class of "topologically similar" 
knots which is, however, much broader than the class represented by the polynomial invariant 
X{t). This enables us to introduce the smoothed measures and distribution functions for rj. 

(ii) The knot complexity rj describes correctly (at least from the physical point of view) 
the limit cases: 77 = corresponds to "weakly entangled" trajectories whereas rj N 
matches the system of "strongly entangled" paths. 

(iii) The knot complexity keeps all nonabelian properties of the polynomial invariants. 

(iv) The polynomial invariant can give exhaustive information about the knot topology. 
However when dealing with statistics of randomly generated knots, we frequently look for 
rougher characteristics of "topologically different" knots. A similar problem arises in statis- 
tical mechanics when passing from the microcanonical ensemble to the Gibbs one: we lose 
some information about details of particular realization of the system but acquire smooth- 
ness of the measure and are able to apply standard thermodynamic methods to the system 
in question. 

The main purpose of the present section is the estimation of the limit probability distri- 
bution of rj for the knots obtained by randomly generated closed -B„,-braids of the length A^. 
It should be emphasized that we essentially simplify the general problem "of knot entropy". 
Namely, we introduce an additional requirement that the knot should be represented by a 
braid from the group i?„ without fail. 

We begin the investigation of the probability properties of algebraic knot invariants by 
analyzing statistics of the random loops ( "Brownian bridges" ) on simplest non- commutative 
groups. Most generally the problem can be formulated as follows. Take the discrete group 
Qn with a fixed finite number of generators {5^1, . . . , Qn-i}- Let u be the uniform distribution 
on the set {gi, . . . , gn-i, gi^, ■ ■ ■ , dn-i}- convenience we suppose hj = gi for j = i and 
hj = g~^ ioT j = i + n — 1] vijij) = for any j. We construct the (right-hand) side 
random walk (the random word) on Qn with a transition measure u, i.e. the Markov chain 
{^n}, ^0 = e G and Prob(^j = = v) = u^v'^u) = 2^^- It means that with the 

probability we add the element haj^ to the given word Hn-i = haiha2 ■ ■ ■ ha^.i from 
the right-hand side 0. 



"^Analogously we can construct the left-hand side Markov chain. 
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The random word W formed by N letters taken independently with the uniform probability 
distribution v = from the set {gi, . . . , g^-i, gi^ , . . . , Qn-i} is called the Brownian 

bridge (BB) of length N on the group Qn if the shortest fprimitivej word ofW is identical 
to the unity. 

Two questions require most of our attention: 

1. What is the probabihty distribution P{N) of the Brownian bridge on the group Qn- 

2. What is the conditional probabihty distribution P{k,m\N) of the fact that the sub- 
word W consisting of first m letters of the A^-letter word W has the primitive path k under 
the condition that the whole word W is the Brownian bridge on the group Qn- (Hereafter 
P{k,m\N) is referred to as the conditional distribution for BB.) 

It has been shown in the paper that for the free group the corresponding problem can 
be mapped on the investigation of the random walks on the simply connected tree. Below we 
represent shortly some results concerning the limit behavior of the conditional probability 
distribution of BB on the Cay ley tree. In the case of braids the more complicated group 
structure does not allow us to apply the same simple geometrical image directly. Nevertheless 
the problem of the limit distribution for the random walks on i?„ can be reduced to the 
consideration of the random walk on some graph C(r). In case of the group we are able 
to construct this graph evidently, whereas for the group i?„ {n > 4) we give upper estimations 
for the limit distribution of the random walks considering the statistics of Markov chains on 
so-called local groups. 



2. Random process on PSL{2,W.), and limit distribution of powers of Alexander invariant 

We begin with computing the distribution function for the conditional random process on 
the simplest nontrivial braid group B^. The group B^ can be represented by 2 x 2 matrices. 



To be specific, the braid generators ai and (T2 in the Magnus representation look as 
follows: 

..^(-;;); (.n.i6) 

where t is "the spectral parameter". It is well known that for t = —1 the matrices (Xi and o"2 
generate the group PSL{2,'Z.) in such a way that the whole group B^ is its central extension 
with the center 

J.6A Q 

(a.a^a,)'' = (aacnas)^" = (aids)^" = (^2^1)'^^ = | ^ J (III.17) 

ij 

First restrict ourselves with the examination of the group PSL{2,'Z), for which we define 
CTi = o"i and a2 = o'2 (at t = —1). 

The canonical representation of PSL{2,'Z.) is given by the unimodular matrices S,T: 
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The braiding relation cricr20'i = (T2<^iO'2 in the {S*, T}-representation takes the form 

S^TS'^T-^ = 1 (III.19) 

in addition we have 

= {STf = 1 (III.20) 

This representation is well known and signifies the fact that in terms of {S, T}-generators 
the group S'L(2,Z) is a free product Z"^ ® of two cyclic groups of the 2nd and the 3rd 
orders correspondingly. 

The connection of {S,T} and {5"i,cr2} is as follows 

a,=T (T = ai) 

(III.21) 

5-2 = T'^ST'^ {S = ox020\) 

The modular group PSL{2,1.) is a discrete subgroup of the group PS'L(2,R). The 
fundamental domain of PSL{2,'Z) has the form of a circular triangle ABC with angles 
|o, |, || situated in the upper half-plane lm( > of the complex plane C = ^ + ^^7 (see fig.|^ 
for details). According to the definition of the fundamental domain, at least one element 
of each orbit of PSL[2,'Z.) lies inside ABC-domain and two elements lie on the same orbit 
if and only if they belong to the boundary of the yli?C-domain. The group P5'L(2,Z) is 
completely defined by its basic substitutions under the action of generators S and T: 

S: C ^ -1/C 

(III.22) 

T: + 1 



fundamental domain 
of the modular group 




FIG. 8. The Riemann surface for the modular group The graph C{T) representing the topological 
structure of PSL{2, 2) is shown by the dashed line. 
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Let us choose an arbitrary element Co from the fundamental domain and construct a 
corresponding orbit. In other words, we raise a graph, C(r), which connects the neighboring 
images of the initial element Co obtained under successive action of the generators from the 
set {S,T, S~^,T~^} to the element Co- The corresponding graph is shown in the fig.|| by 
the broken line and its topological structure is clearly reproduced in fig.^. It can be seen 
that although the graph C(T) does not correspond to the free group and has local cycles, 
its "backbone", C(7), has Cayley tree structure but with the reduced number of branches 
as compared to the free group C(T2). 




FIG. 9. The graph C{T) and its backbone graph C{j) (see the explanations in the text). 

Turn to the problem of limit distribution of a random walk on the graph C(r). The walk 
is determined as follows: 

1. Take an initial point ("root") of the random walk on the graph C(r). Consider 
the discrete random jumps over the neighboring vertices of the graph with the transition 
probabilities induced by the uniform distribution u on the set of generators {ai, a2, 5"f ^, 5"^^}. 
These probabilities are (see Eq.( [lII.21|) ) 

Prob(a = TCo I ^n-i = Co) = i 
Prob(a = (T-iST-i)Co I a-i = Co) = i 

(111.23) 

Prob(a = T-^Co I ^n-i = Co) = i 
Prob(a = {TS-'T)Co I Cn-i = Co) = i 

The following facts should be taken into account: the elements S(q and S'~^Co represent one 
and the same point, i.e. coincide (as it follows from Eg. ( |111.2"2D ); the process is Markovian 
in terms of the alphabet {5"i, . . . , 5"^^} only; the total transition probability is conserved. 

2. Define the shortest distance, k, along the graph between the root and terminal points 
of the random walk. According to its construction, this distance coincides with the length 
|W^{s,T}| of the minimal irreducible word W{s,t} written in the alphabet {S,T, ,T~^}. 
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The link of the distance, k, with the length | ^{^1,5-2} | of the minimal irreducible word ^^{5-1,52} 
written in terms of the alphabet {a"i, 5"2, 5"f ^, a^^} is as follows: (a) |W^{5i,5-2}| = if and 
only if k = 0; (b) for A; 1 the length |Vr{5-^_s2}l has asymptotic: |Vr{5-^^CT2}| = k + o{k). 
We define the "coordinates" of the graph vertices in the following way (see fig.P): 

(a) We apply the arrows to the bonds of the graph F corresponding to T-generators. The 
step towards (backwards) the arrow means the application of T (T~^). 

(b) We characterize each elementary cell of the graph F by its distance, /x, along the 
graph backbone 7 from the root cell. 

(c) We introduce the variable a = {1,2} which numerates the vertices in each cell only. 
We assume that the walker stays in the cell M located at the distance /i along the backbone 
from the origin if and only if it visits one of two in-going vertices of M. Such labelling gives 
unique coding of the whole graph C(F). 

Define the probability Ua{fi, N) of the fact that the A^-step random walk along the graph 
C(F) starting from the root point is ends in a- vertex of the cell on the distance of /i steps 
along the backbone. It should be emphasized that f/Q,(/i, N) is the probabihty to stay in any 
of Afry{fi) = 3 ■ 2^^^ cells situated at the distance ^ along the backbone. 

It is possible to write the closed system of recursion relations for the functions f/a(/i, A^). 
However, here we attend to rougher characteristics of random walk. Namely, we calculate 
the "integral" probability distribution of the fact that the trajectory of the random walk 
starting from an arbitrary vertex of the root cell O has ended in an arbitrary vertex point of 
the cell M situated on the distance /i along the graph backbone. This probability, U{fi, N), 
reads 

uif,,N) = ^ y: t/.(/i,iv) 

^ a={l,2} 

The relation between the distances k, along the graph F, and /i along its backbone 7 is such: 
k = 11 + o{fi) for /i ^ 1, what ultimately follows from the constructions of the graphs C(F) 
and C(7). 

Suppose the walker stays in the vertex a of the cell M located at the distance ^ > 1 from 
the origin along the graph backbone C{'y). The change in /i after making of one arbitrary 
step from the set {ai, 5"2, af ^, a^^} is summarized in the following table: 



a = 1 


a = 2 


a, = T 


|U ^ /i + 1 


ai = T 




d2 = T-^ST-^ 




02 = T-^ST-^ 


fi ^ fi + 1 




fi —>■ fi — 1 


af' = 


^ fi + 1 


a^^ = TS-^T 


|U — > /i + 1 


^2"^ = TS-^T 





It is clear that for any value of a two steps increase the length of the backbone, /i, one 
step decreases it and one step leaves fi without changes. 

Let us introduce the effective probabihties: pi - to jump to some specific cell among 
3 neighboring ones of the graph C(F) and p2 - to stay in the given cell. Because of the 
symmetry of the graph, the conservation law has to be written as 3pi +p2 = 1. By definition 
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we have: pi 



V 



\. Thus we can write the following set of recursion relations for the 



integral probability ?7(/i, A^) 

f/(/x, iV + 1) = if/(/i + 1, iV) + if/(/i, iV) + if/(/i - 1, N) (/i > 2) 
\J{ix, iV + 1) = if/(/i + 1, iV) + if/(/i, N) ifi = 1) 

[/(/i,iV = 0) = 5^,i 



^111.24) 



The solution of Eq.(|TL2^) gives the limit distribution for the random walk on the group 
PLS{2,Z). 

The probability distribution U {k, N) of the fact that the randomly generated N -letter word 
^{5-i,5-2} '"'^^^ uniform distribution y = \ over the generators {5"i, 5-2, af ""^j 0'2^^} can be 
contracted to the minimal irreducible word of length k, has the following limit behavior 



U{k,N) 



h 



h 



N 



1 



^/^{A-h) \A{h-2) 



iV3/2 
k 



2^/2 



exp 



k = 
1 < k 



(III.25) 



where h 



2 + ^. 

2 



Corollary 1 The probability distribution U{k, m\N) of the fact that in the randomly gener- 
ated N -letter trivial word in the alphabet {5"i, o"2, (T2 ^} the sub-word of first m letters 
has a minimal irreducible length k reads 



U{k,m\N) 



h 



k^ 



v/¥(4 - h) {m{N - m))3/2 



exp 



k^h 



m N — m 



(III.26) 



Actually, the conditional probability distribution U{fi,m\N) that the random walk on 
the backbone graph, C{'j), starting in the origin, visits after first m = const) steps some 
graph vertex situated at the distance /i and after steps returns to the origin, is determined 
as follows 



where TV^ = 3 ■ 2^-^) and U{n, N) is given by ([nT25| ). 



The problem considered above helps us in calculating the conditional distribution func- 
tion for the powers of Alexander polynomial invariants of knots produced by randomly 
generated closed braids from the group B3. 

Generally the closure of an arbitrary braid b G -B„ of the total length gives the knot 
(link) K. Split the braid b in two parts b' and b" with the corresponding lengths m and 
N — m and make the "phantom closure" of the sub-braids b' and b" as it is shown in fig.pil|. 
The phantomly closed sub-braids b' and b" correspond to the set of phantomly closed parts 
( "sub- knots" ) of the knot (link) K. The next question is what the conditional probability 
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to find these sub-knots in the state characterized by the complexity t] when the knot (hnk) 
as a whole is characterized by the complexity rj = (i.e. the topological state of K "is 
close to trivial"). 





FIG. 10. Construction of Brownian bridge for knots. 
Returning the the group B^, introduce normalized generators as follows 

(det 



\<J 



±ii 



To neglect the insignificant commutative factor dealing with norm of matrices ai and a2- 
Now we can rewrite the power of Alexander invariant (Eq.( p.ll.I^ )) in the form 



(III.28) 



where #(+) and ) are numbers of generators cTq,^ or a^^ in a given braid and t] is the 
power of the normalized matrix product Ylf=i WcajW- The condition of Brownian bridge 
implies 77 = (i.e. #(+) — #(— ) = and rj = 0). 
Write 



llaill =r(t); \\a2\\=T-\t)S{t)T-\t) 
where T(t) and S(t) are the generators of the "t-deformed" group PSLt{2, 



(III.29) 



m 

S{t) 









(-t)-V2 

(-t)i/2 



1 (-t)-l 

1 

1 
-1 



(III.30) 



The group PSLt{2,'Z) preserves the relations of the group P5'L(2,Z) unchanged, i.e., 
iT{t)S{t)f = S\t) = T{t)S^{t)T-\t)S~^it) = 1 (compare to Eq.(|ITLT9D). Hence, if we 
construct the graph C{rt) for the group PSLt{2,'Z) connecting the neighboring images of an 



39 



arbitrary element from the fundamental domain, we ultimately come to the conclusion that 
the graphs C(Tt) and C(T) (fig.^ are topologically equivalent. This is the direct consequence 
of the fact that group is the central extension of PSL{2,'Z.). It should be emphasized that 
the metric properties of the graphs C{rt) and C(r) differ because of different embeddings 
of groups PSLt{2,'Z) and PSL{2,'Z.) into the complex plane. 

Thus, the matrix product Tl^Li ll'^ojll for the uniform distribution of braid generators 
is in one-to-one correspondence with the iV-step random walk along the graph C(r). Its 
power coincides with the respective geodesies length along the backbone graph C{'-f). Thus 
we conclude that limit distribution of random walks on the group -B3 in terms of normalized 
generators ( [111.29| ) is given by Eq.( [111.2"5| ) where k should be regarded as the power of the 
product na=i ||<7aj||- Hence we come to the following statement. 

Take a set of knots obtained by closure of B^-braids of length N with the uniform dis- 
tribution over the generators. The conditional probability distribution U{rj,m\N) for the 
normalized complexity T] of Alexander polynomial invariant (see ( 111.28| )) has the Gaussian 
behavior and is given by Eq. (pi.26| ) where k = rj. 



B. Random walks on locally free groups 

We aim at getting the asymptotic of conditional limit distributions of BB on the braid 
group Bn- For the case n > 3 it presents a problem which is unsolved yet. However we can 
estimate limit probability distributions of BB on Bn considering the limit distributions of 
random walks on the so-called "local groups" ( 



The group CJ^n+i{d) we call the locally free if the generators, {/i, . . . , /„} obey the 
following commutation relations: 

(a) Each pair {fj, fk) generates the free subgroup of the group JF„ if |j — /c| < 

(b) fjfk = fkfj for \j -k\>d. 

(Below we restrict ourselves to the case d = 2 where CJ-'n+i{2) = CJ-'n+i)- 
The limit probability distribution for the N-step random walk {N ^ 1) on the group Tn+i 
to have the minimal irreducible length fi is 



;iH.31) 



We propose two independent approaches valid in two different cases: (1) for n = 3 and 
(2) for n > 1. 

(1) The following geometrical image seems useful. Establish the one-to-one correspon- 
dence between the random walk in some n-dimensional Hilbert space £?i"(xi, . . . and 
the random walk on the group CJ^n+i, written in terms of generators {/i, . . . , /~^}. To be 
more specific, suppose that when a generator, say, fj, (or f~^) is added to the given word 
in CJ-'n, the walker makes one unit step towards (backwards for f^^) the axis [0,Xj[ in the 
space CTC^Xi, . . . , 

Now the relations (a)-(b) of the definition of the locally free group could be reformulated 
in terms of metric properties of the space £7i". Actually, the relation (b) indicates that 
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successive steps along the axes [0,Xj[ and [0,Xfc[ (|j — k\ > 2) commute, hence the section 
{xj,Xk) of the space CTC"' is flat and has the Euclidean metric dx'j + dx\. Situation with 
the random trajectories in the sections {xj,Xj±i) of the Hilbert space CT-C appears to be 
completely different. Here the steps of the walk obey the free group relations (a) and the 
walk itself is mapped to the walk on the Cayley tree. It is well known that Cayley tree can 
be uniformly embedded (without gaps and selfintersections) into the 3-pseudosphere which 
gives the representation of the non-Euclidean plane with the constant negative curvature. 
Thus, sections {xj,Xj+i) have the metric of Lobachevskii plane which can be written in the 
form ^{dxj + dx'j_^_^). 

For the group CJ-'^ these arguments result in the following metric of appropriate space 

ds' = + + dxl 

xi 

Actually, the space section {xi,X3) is flat whereas the space sections {xi,X2) and {x2,X3) 
have Lobachevskii plane metric. The noneuclidean (hyperbolic) distance between two points 
M' and M" in the space is defined as follows 

coshf^iM'M") = 1 + E (^.(M') - x.{M")f (111.33) 

where {xi,X2, x^} are the euclidean coordinates in the 3D-halfspace X2 > and ^ is regarded 
as geodesies on a 4-pseudosphere (Lobachevskii space). 

Some well known results concerning the limit behavior of random walks in spaces of 
constant negative curvature are reviewed in the next Section where solutions of the diffusion 
equations in the Lobachevskii plane and space are given by Eg. ( [111. 5"^ ) and Eq.( |111.6"0D 
correspondingly. Thus we can conclude that the distribution function for random walk in 
Lobachevskii space Vs{fi,N) defined by Eqs. ([1II.60| )-( |1II.63D gives also the probability for 
the iV-letter random word (written in terms of uniformly distributed generators on JF4) to 
have the primitive word of length /x (see Eg. ([III. 3 1|) ). 

(2) For the group LTn+\ {n ^ 1) we extract the limit behavior of the distribution 
function evaluating the probabilities to increase and to decrease the length of the primitive 
word if we randomly add one extra letter to the given word. We follow below the line 
proposed by J. Desbois [0,0] 



Let us point out the main steps of our computations: 

(a) We generate randomly (with uniform probability distribution) the words of lengths 

G [1000; 20000], while the number of generators, n, varies in the interval [3; 200]. 
The number of randomly generated words is of order of 1000. 

(b) We reduce the given word till the minimal irreducible (primitive) word. This can be 
done by using the braid (or locally free) group relations. The numerical procedure is 
as follows. First, we try to push each braid generator in the word as far as possible to 
the left. Some reductions can occur after that. Then, we play the same game but in 
the opposite direction, pushing each braid generator to the right performing possible 
reductions of the word, then — to the left again and so on... If no reductions occur 
during two consecutive steps, we stop the process. 



41 



We compute the following quantities for braid and locally free groups: 
The mean length of the shortest (primitive) word (/i) 



oo 

(/i) = '-^ (111.34) 

E Z{fi,N) 

/i=0 

and the variance Var(/i) 

oo 

Var(/i) ^ (/i^) - = tzl (111.35) 

E ^(/i,iV) 

fj.=0 

The results of numerical simulations for the word statistics on braid (-Bn) and locally 
free {CJ^n{d)) groups are presented in the Table 1. 

Table iQ. 



Groups 


Bn 










N 


Var{/i) 
N 


N 


Var{/i) 
N 


N 


Var(At) 
N 


N 


Var(At) 
N 


n=3 


0.29 


0.85 


0.50 


0.76 


0.50 


0.76 


0.50 


0.75 


n=5 


0.49 


0.77 


0.60 


0.63 


0.71 


0.48 


0.75 


0.46 


n=10 


0.56 


0.63 


0.65 


0.56 


0.77 


0.40 


0.82 


0.34 


n=20 


0.59 


0.63 


0.66 


0.54 


0.79 


0.39 


0.84 


0.29 


n=50 


0.61 


0.61 


0.67 


0.56 


0.80 


0.38 


0.85 


0.27 


n=100 


0.61 


0.61 


0.67 


0.52 


0.80 


0.36 


0.86 


0.26 


11=200 


0.61 


0.60 


0.67 


0.53 


0.80 


0.35 


0.86 


0.26 



The maximal standard deviations in the Table 1 (and everywhere below) are: 

±0.01 for the mean value (/i) /N 
±0.05 for the variance Var(/i)/A^ 

C. Analytic Results for Random Walks on Locally Free Group 

Let us estimate now the quantities (//) /N and Yai{^)/N analytically. We present below 
two different approaches called "dynamical" and "statistical". The "dynamical" approach 
is based on simple estimation of the probability to reduce the primitive word by random 
adding one extra letter. The estimate obtained by this method is in very good agreement 



^The groups CJ-n{d) are completely free when d>n — 1. 
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with corresponding numerical simulations. However the "statistical" approach dealing with 
rigorous enumeration of all nonequivalent primitive words in the locally free group CJ^n{d) 
leads to another answer. The rest of this section is devoted to the explanation of the 
abovementioned discrepancy. 

Dynamical Consideration. Under the conditions 

^ \ (111.36) 

we can easily develop the dynamical arguments which are in rather good agreement with 
the results of numerical simulations presented above. The last inequality in ( |1II.36|) ensures 



the conditions, sufficient for finding the limit probability distribution of Markov chains on 
the groups of n generators. Actually, the number of letters in the word, should be much 
larger that the number of all possible pairs in the set of 2n letters^. Only in this case the 
corresponding Markov process has the reliable distribution function. The number of pairs 
is of order 4?t,^, so we arrived at the inequality stated in ( [ill.36| ). 



Take a randomly generated A^-letter word W. This word is characterized by the length 
of the primitive word Wp (recall that Wp is the length of the word W obtained after all 
possible contractions allowed by the structure of the group 

Let us compute the probability n{d) of the fact that the primitive word Wp will be 
shortened in one letter after adding of the letter fi {i G [1,''^]) to the word W from the 
right-hand side. It is easy to understand that the primitive word Wp can be reduced if: 



a) The last letter in the word Wp is just ^. The probability of such event is ^ 



2n' 



b) The letter before the last in the word Wp is ^ and the last letter commutes with the 
letter /j. The probability of such event is ^ (^1 ^'^"^ 



2n 

c) The third letter from the right end of the word Wp is and two last letters commute 
with the letter /j. The probability of such event is 2^ (l — 

d) ... and so on. 

Finally we arrive at the following expression for the probability TT{d): 

M<i)^ig(l-V) (111.37) 

The procedure described above assumes that the letters remaining in the word Wp are 
uniformly distributed — as in the initial (nonreduced word W). The absence of "boundary 
effects" is ensured by the condition ( [II.36|) . 



^The total number of generators is 2n because each of n generators has the inverse one. 
^Our consideration is valid for any values of d. 
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Once having the probabihty vr((i), we can write down the master equation for the proba- 
bihty P{fi, N) of the fact that in randomly generated A^-letter word the primitive path has 
the length fj, 

P{lj,N + l) = {l-7r{d))P{i2-l,N)+7r{d)P{i2 + l,N) (/i > 2) (111.38) 

where the relation between N) and the partition function Z{fi, N) introduced above is 
as follows 

fj.=0 

The recursion relation ( |111.38| ) coincides with the equation describing the random walk on 
the halfline with the drift from the origin or, what is the same, with the equation describing 
the random walk on the simply Cayley tree with the coordinational number 

2eff = =Ad-2 (III.39) 

Taking into account the last analogy we can complete the Eq. (|lII.38| ) by the boundary 
conditions 

P(/i = 1, AT + 1) = P(/i = 0, N) + Tx{d) P(/i = 2, N) 

P{lJ = 0,N +1) = nP{fi = l,N) (III.40) 
P(/i,iV = 0)=5^,o 

It is noteworthy that these equations are written just for the Cayley tree with z^s branches. 
The actual structure of the graph corresponding to the group CJ-'n{d) is much more com- 
plex, thus Eqs.( [lll.l0D should be regarded as an aproximation. However the exact form of 
boundary conditions does not influence the asymptotic solution of Eq. (|III.38|) in vicinity of 
the maximum of the distribution function: 



P(^, N) ^ ^ exp { , "^"^ , ( fi - ^^^iv) ' \ (III.41) 

2^27r(zefr - l)N \ 8(zeff-l)iV^^ . ) ( ^ > 



Thus, we find 



^ Zeff - 2 _ 2rf- 2 
Var(/i, d) A{zcs - 1) M 



[IIIA2) 



N zl^ {2d -If 

Substituting in Eq.( [111.42D (i = 2, 3, 4 we get the following numerical values: 

for = 2 

for (i = 3 

for = 4 





2 


Var(/i, d) 


5 


N 


~ 3' 


N 


~ 9 




4 


Var(/i, d) 


9 


N 


~ 5' 


N 


~ 25 


W)) 


6 


Var(/i, d) 


13 


N 


~ 7' 


N 


~ 49 
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what is in the excellent agreement with the asymptotic values (n ^ 1) from the Table 1 for 
the same groups. 

Another statistical problem appears when we are interested in the consideration of the 
target space of the group LTn+i-, i-e. in the evaluation of the number of nonequivalent 
primitive words in the froup CJ-'n+i (see for details ||53[]). 

Let Vn{^i) be the number of all nonequivalent primitive words of length /i on the group 
CJ-'n+i- When /i 1, has the following asymptotic: 



= const 



2 3- 



471 



2' 



'n? 



n > 1 (III.43) 



To get Eq. ( [111.4^ ) we write each primitive word Wp of length /i in the group LTn+x in 
the so-called normal order (all /q. are different) similar to so-called "symbolic dynamics" 
used in consideration of chaotic systems 

w^p = (/..r (/..r---(/«.r (ni.44) 



where ^ |mj| = [i [rrii ^ Q'i i; I < s < fi) and sequence of generators /„. in Eq.( p.ll.4^ ) ft 



or 



i=l 



all satisfies the following local rules: 

(i) If fa, = /i, then e {/2, /s, • • • /n-i}; 

(ii) If fa, = fk (l < k <n - 1), then fa,+^ e {fk-i, fk+i, ■ ■ ■ fn-i}; 

(iii) If fa, = fn, then fa,+^ = fn-l- 

These local rules prescribe the enumeration of all distinct primitive words. If the sequence 
of generators in the primitive word Wp does not satisfy the rules (i)-(iii), we commute 
the generators in the word Wp up the normal order is restored. Hence, the normal order 
representation provides us with the unique coding of all nonequivalent primitive words in 
the group CJ-'n+i- 

The calculation of the number of distinct primitive words, V^(/i), of the given length ^ 
is rather straightforward: 



K(/i) = E^(^) E' ^ 



{mi,...,ms} Li=l 



s 



(III.45) 



where R{s) is the number of all distinct sequences of s generators taken from the set 
{/i,...,/„} and satisfying the local rules (i)-(iii) while the second sum gives the num- 
ber of all possible representations of the primitive path of length /i for the fixed sequence 
of generators ("prime" means that the sum runs over all rrii ^ for 1 < i < s; A is the 
Kronecker A-function). 

It should be mentioned that the local rules (i)-(iii) define the generalized Markov chain 
with the states given by the n x n coincidence matrix T„ where the rows and columns 
correspond to the generators /i, . . . , /„: 



45 



fn{d) 





fl 


/2 


/3 


/4 




fn-l 


fn 


h 





1 


1 


1 




1 


1 




1 





1 


1 




1 


1 


fs 


1 


1 





1 




1 


1 


u 





1 


1 







1 


1 


















fn-l 



















1 


fn 
















1 
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The number of all distinct normally ordered sequences of words of length s with 
allowed commutation relations is given by the following partition function 



Rn{s,d) = Vin fn{d) 



'out 



(III.47) 



where 











> 


n 






1 




Vi„ = (rTT^~T) 


and 


Vout = 


1 


> 













(III.48) 



Supposing that the main contribution in Eq.( |111.45D results from s 1 we take for Rn{s) 
the following asymptotic expression 



Rnis] 



(a: 



max\s 



47r2 ^ / 1 
3-^ + 



(III.49) 



where A^"^ is the highest eigenvalue of the matrix T„ (n ^ 1). 

The remaining sum in Eq.( |111.4"5D is independent of -R(s), so its calculation is trivial: 



{mi,...,r?is} 



i=l 



l)!(/i 



(III.50) 



Collecting all terms in Eq. ([lII.45| ) and evaluating the sum over s we arrive at Eq. (pi.43| ). 
The value c?) is growing exponentially fast with and the "speed" of this grows is 

clearly represented by the fraction 



7 



^111.51) 



where z^^ is the coordinational number of effective tree associated with the locally free group. 

Thus, the random walk on the group CJ^n+i can be viewed as follows. Take the free 
group r„ with generators {/i, . . . , /„} where all fi {1 < i < n) do not commute. The group 
r„ has a structure of 2n-branching Cayley tree, C(r„), where the number of distinct words 
of length /i is equal to Vn^jj), 
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VM = 2n{2n - 1)^"^ (III.52) 

The graph C{CJ-'n+i) corresponding to the group CJ-'n+i can be constructed from the graph 
C(r„) in accordance with the following recursion procedure: (a) Take the root vertex of the 
graph C(r„) and consider all vertices on the distance fi = 2. Identify those vertices which 
correspond to the equivalent words in group CJ-'n+i] (b) Repeat this procedure taking all 
vertices at the distance /i = (1,2,.. .) and "gluing" them at the distance /i + 2 according 
to the definition of the locally free group. By means of the described procedure we raise 
a graph which in average has Zes — 1 distinct branches leading to the "next coordinational 
sphere" . Thus this graph coincides (in average) with Zeff-branching Cayley tree. 

Although the local structure of the graph C (CJ-'n+i) is very complex, Eq. (|lII.51| ) enables 
us to find the asymptotic of the random walk on the graph C (CJ-'n+i)- Once having Zes, we 
can write down the master equation for the probability V{fi, N) to find the walker at the 
distance fi from the origin after random steps on the graph C (CJ-'n+i) 

V{fi,N + l) = {l-—)P{fi-l,N) + —P{fi + l,N) (/i>2) (III.53) 



The recursion relation (|III.53| ) coincides with the equation describing the random walk on the 
half-line with the drift from the origin. Taking into account this analogy we can complete the 
Eq.( [111.5"3| ) by the boundary conditions ||5^. However the exact form of boundary conditions 



does not influence the asymptotic solution of Eq.(|TL5^) in vicinity of the maximum of the 
distribution function: 

N) ^ ^ exp I -^Sl^ ( f, - M^iV 



2J27i{zes - l)N I 8(^cfr - 1)A^ V z^s 



Thus we obtain the desired distribution function (Eq.( |III.3TD ) for the primitive word length 
for the random walk on the group CTn+\- 

The Eq.( |111.3"l| ) gives the estimation from below for the limit distribution of the primitive 
words on the group Bn for n 3> 1 . 

We find further investigation of the random walks on the groups CJ^n+i{d) for different 
values of d very perspective. It should give insight for consideration of random walk statistics 
on "partially commutative groups" . Moreover, the set of problems considered there has deep 
relation with the spectral theory of random matrices. 



D. Brownian bridges on Lobachevskii plane and products of non-commutative 

random matrices 

The problem of word enumeration on locally non-commutative group has evident connec- 
tion with the statistics of Markov chains on graphs having the Cayley tree-like structure and, 
hence, with random walk statistics on the surfaces of a constant negative curvature. (We 
stressed once that the Cayley tree-like graphs are isometrically embedded in the surfaces of 
a constant negative curvature). 

Recall that the distribution function, P(r,t), for the free random walk in D-dimensional 
Euclidean space obeys the standard heat equation: 

|-P(r,t)=PAP(r,t) 
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with the diffusion coefficient V = and appropriate initial and normahzation conditions 



P(r,t = 0) = 5(r) 
j P{r,t)dr = 1 

Correspondingly, the diffusion equation for the scalar density P(q, t) of the free random 
walk on a Riemann manifold reads (see [B^ for instance) 



t) 



(111.54) 



where 



P(q,t = 0) = 5(q) 
^P{q,t)dq=l 



(III.55) 



and Qik is the metric tensor of the manifold; g = det Qik- 

Eq. (|III.54|) has been subjected to thorough analysis for the manifolds of the constant 
negative curvature. Below we reproduce the corresponding solutions for the best known 
cases: for 2D- and 3D-Lobachevskii spaces (often referred to as 3- and 4-pseudospheres) 
labelling them by indices "p " and "s " for 2D- and 3D-cases correspondingly. 

For the Lobachevskii plane one has 



\9ik\ 



1 








sinh^ fx 



(III.56) 



where fi stands for the geodesies length on 3-pseudosphere. The corresponding diffusion 
equation now reads 



d 

—Pp{fi,^,t) = V 



, a 1 

+ coth fi— — h 



Pp{fX,(fi,t) 



^ dfi^ ' ^ dfi ' sinh^ fi dip^ j 
The solution of Eq.(|TL5^) is believed to have the following form 

eexp(-^) 



(III.57) 



— *2 

e 4 



a/ cosh ^ — cosh /i 



d^ 



e 4 



1/2 



(III.58) 



exp 



fi^ 



AW 



ATrtV \ sinh fi ^ 

For the Lobachevskii space the corresponding metric tensor is 

10 
\9ik\ \ = sinh^/x 

sinh^ fi sin^ 9 



[111.59) 
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Substituting Eq. (|111.54|) for Eq. (pL59| ) we have 



P,(/i,t) = —^^^=-J^exp (-^] (III.60) 



87r^7r(tP)3sinh/i V AW 
For the first time this spherically symmetric solution of the heat equation (Eq.( [111.5^ ) in 



the Lobachevskii space was received in |63 



In our opinion one fact must be given our attention. The distribution functions Pp(/i,t) 
and Ps(/^, t) give the probabilities to find the random walk starting at the point /i = after 
time t in some specific point located at the distance /x in corresponding noneuclidean space. 
The probability to find the terminal point of a random walk after time t somewhere at the 
distance is 

Pp,,(/i, t) = Pp,s(/i, t)Af^M (III.61) 

where 

A/'p(/i) = sinh/i (III.62) 
is the perimeter of circle of radius fi on the Lobachevskii plane and 

A/;(/i) = sinhV (III.63) 

is the area of sphere of radius fi in the Lobachevskii space. 

The difference between Pp and Vp^s is insignificant in euclidean geometry, whereas in 
the noneuclidean space it becomes dramatic because of the consequences of the behavior of 
Brownian bridges in spaces on constant negative curvature. 

Using the definition of the Brownian bridge, let us calculate the probabilities to find 
the A^-step random walk (starting at fi = 0) after first t steps at the distance /i in the 
Lobachevskii plane (space) under the condition that it returns to the origin on the last step. 
These probabilities are [N — > oo) 

Vp{^,t\0,N) = = — -fiexpi-— - + 



Pp(0,t) ATTVt{N-ty [ AV\t N-t 



P,{0,t) 87rt3/2(Ar-t)3/2^ ^\ 4V\t N-t 

(III.64) 

Hence we come to the standard Gaussian distribution function with zero mean. 

Equations ( [1II.64|) describing the random walk on the Riemann surface of constant nega- 



tive curvature have direct application to the conditional distributions of Lyapunov exponents 
for products of some non-commutative matrices. Let us consider the first of Eqs. ( [111.6^ ). 



Changing the variables /x = In Y^jfi ! V = arg^; where z = x + iy; z = x — iy we map the 
3-pseudosphere (/i, ip) onto the unit disk < 1 known as the Poincare representation of the 
Lobachevskii plane. The corresponding conformal metric reads dP = ^j-^^^jT- Using the con- 

formal transform z = we recover the so-called Klein representation of Lobachevskii 

1 — iw 

plane, where dl'^ = and the model is defined in Imw > (w = u + iv: w = u — iv). 
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The following relations can be verified using conformal representations of the Lobachev- 
skii plane metric (see, for instance, fl^])- The fractional group of motions of Lobachevskii 
plane is isomorphic to: 

(i) the group SU{1, 1)/±1 = PSU{1, 1) in the Poincare model; (ii) the group SL{2, R)/± 
1 = PSL{2,R) in the Klein model. 



Moreover, it is known (see, for example, |0) that the Lobachevskii plane H can be iden- 



tified with the group SL{2,R)/ S0{2). This relation enables us to resolve (at least qualita- 
tively) the following problem. Take the Brownian bridge on the group 7i = SL{2, R) / 5*0(2), 
i.e. demand the products of N independent random matrices A^fcG7i(0</c< N) to 
be identical to the unit matrix. Consider the limit distribution of the Lyapunov exponent, 
6, for the first m matrices in that products. To have a direct mapping of this problem on 
the random walk in the Lobachevskii plane, write the corresponding stochastic recursion 

equation for some vector = [ ) 

\VkJ 

Wfc+i = MkW,; Wo = Q (in.65) 

where A^^ G 7i for all k G [1, A^]- The BB-condition means that 

Wn = Wi for iV > 1 (IIL66) 

Let us consider the simplest case 

Mk = l + Mk, norm[Mfc] < 1 (IIL67) 

In this case the discrete dynamic equation ( [1I1.65| ) can be replaced by the differential one. 
Its stationary measure is determined by the corresponding Fokker-Plank equation ([111.541) . 



The Lyapunov exponent, S of product of random matrices A4 coincides with the length of 
geodesies in the Klein representation of the Lobachevskii plane. Hence, under the conditions 
( [111. 66] ), ( [111.67[ ) we have for 6 the usual Gaussian distribution coinciding with the first of 



Eq. ( [111.6^ ). Without the BB-condition (i.e. for "open walks") we reproduce the standard 
Fiirstenberg behavior [^]. 

Although this consideration seems rather crude (for details see Appendix A), it clearly 
shows the origin of the main result: 

The "Brownian bridge" condition for random walks in space of constant negative curva- 
ture makes the space "effectively flat" turning the corresponding limit probability distribution 
for random walks to the ordinary central limit distribution. 

The question whether this result is valid for the case of the random walk in noneuclidean 
spaces of non-constant negative curvature still remains. 

Finally we would like to introduce some conjectures which naturally generalize our con- 
sideration. 

The complexity rj of any known algebraic invariants (Alexander, Jones, HOMFLY) for the 
knot represented by the Bn-braid of length N with the uniform distribution over generators 
has the following limit behavior: 

^(^' ^) ~ ( -«(^)^ + Pin)v - j^j (111.68) 
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where a{n), (3{n), S{n) are numerical constants depending on n only. 

The knot complexity rj in ensemble of Brownian bridges from the group shown in 
fig.0 has Gaussian distribution, where 

iv) = 0; ) = l5{n)N (III.69) 

These conjectures are to be proven yet. The main idea is to employ the relation between 
the knot complexity 77, the length of the shortest noncontractible word and the length of 
geodesies on some hyperbolic manifold. 



IV. CONFORMAL METHODS IN STATISTICS OF RANDOM WALKS WITH 

TOPOLOGICAL CONSTRAINTS 

The last few years have been marked by considerable progress in understanding the 
relationship between Chern-Simons topological field theory, construction of algebraic knot 
and link invariants and conformal field theory (see, for review, 0]). 



Although the general concepts have been well elaborated in the field-theoretic context, 
their application in the related areas of mathematics and physics, such as, for instance, 
probability theory and statistical physics of chain-like objects is highly limited. 

The present Section is mainly concerned with the conformal methods in statistical anal- 
ysis which allow us to correlate problems discussed in Chapters 1 and 2 and the limit distri- 
butions of random walks on multiconnected Riemann surfaces. To be more specific, we show 
on the level of differential equations how simple geometrical methods can be applied to con- 
struction of non-commutative topological invariants. The latter might serve as nonabelian 
generalizations of the Gauss linking numbers for the random walks on multi-punctured Rie- 
mann surfaces. We also study the connection between the topological properties of random 
walks on the double punctured plane and behavior of four-point correlation functions in 
the conformal theory with central charge c = —2. The developed approach is applied to 
the investigation of statistics of 2D-random walks with multiple topological constraints. For 
instance, the methods presented here allow us to extract nontrivial critical exponents for the 
contractible (i.e., unentangled) random walks in the regular lattices of obstacles. Some of 
our findings support conjectures of Sections 2 and 3 and have direct application in statistics 
of strongly entangled polymer chains (see Section 5). 



A. Construction of nonabelian connections for r2 and PS'L(2,Z) from conformal 

methods 

We analyze the random walk of length L with the effective elementary step a {a = 1) on 
the complex plane z = x + iy with two points removed. Suppose the coordinates of these 
points being Mi {zi = (0,0)) and M2 {z2 = (c, 0)) (c = 1). Such choice does not indicate 
the loss of generality because by means of simultaneous rescaling of the effective step, a, of 
the random walk and of the distance, c, between the removed points we can always obtain 
of any arbitrary values of a and c. 
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Consider the closed paths on z and attribute the generators g\, §2 of some group G to 
the turns around the points Mi and M2 if we move along the path in the clockwise direction 
(we apply (yff ^, (yf^^ for counter-clockwise move) — see fig.llll 



—o 




FIG. 11. (a) — the double punctured complex plane z with two basis loops C\ and C2 enclosing 
points Ml and M2; (b) — the universal covering with fundamental domain corresponding to free 
group V^- The contours Pi and P2 are the images of the loops C\ and C2. 

The question is: what is the probability V) for the random walk of length L on the 
plane z to form a closed loop with the shortest noncontractible word written in terms of 
generators {^gi-, gi-, ■, 9^^^ to have the length /i (see also Chapter 2). 

Let the distribution function P(/i, V) be formally written as a path integral with a Wiener 
measure 
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JO \ ds 



W{gi,g2,gi ,92 \z} - /i 



(IV.l) 



where Z = J P{fi, L)dfi and W{. . .\z} is the length of the shortest word on G as a functional 
of the path on the complex plane. 

Conformal methods enable us construct the connection and the topological invariant W 
for the given group as well as to rewrite Eq. ( |1V.1| ) in a closed analytic form which is solvable 
at least in the limit L ^ 00. 

Let C,{z) be the conformal mapping of the double punctured plane z = x + iy on the 
universal covering C, = ^ + i\. The Riemann surface C, is constructed in the following way. 
Make three cuts on the complex plane z between the points Mi and M2, between M2 and 
(00) and between {00) and Mi along the line Ivnz = 0. These cuts separate the upper 
(Im^ > 0) and lower [Imz < 0) half-planes of z. Now perform the conformal transform of 
the half-plane Imz > to the fundamental domain of the group G{gi,g2} — the curvilinear 
triangle lying in the half-plane lm( > of the plane (. Each fundamental domain represents 
the Riemann sheet corresponding to the fibre bundle above z. The whole covering space ( 
is the unification of all such Riemann sheets. 

The coordinates of initial and final points of any trajectory on universal covering ( de- 
termine ( ^6[)- ('^) -^^^ coordinates of corresponding points on z; (b) The homotopy class 



of any path on z. In particular, the contours on ( are closed if and only if W{gi, (72 1 -2} = 1? 
i.e. they belong to the trivial homotopy class. 

Coordinates of ends of the trajectory on universal covering ( can be used as the topo- 
logical invariant for the path on double punctured plane z with respect to the action of the 
group G. 

Thus, we characterize the topological invariant, Inv(C), of some closed directed path C 
starting and ending in an arbitrary point zq 7^ {zi, Z2,oo} on the plane z by the coordinates 
of the initial, Cin(-2o), and final, Cfin(-2o), points of the corresponding contour P in the covering 
space The contour P connects the images of the point zq on the different Riemann sheets. 
Write InV(^)(C) as a full derivative along the contour C: 

Inv(,) (C) Cin - Cfin = £ ^dz (IV.2) 

dC^ (^) 

The physical interpretation of the derivative — - — is very straightforward. Actually, the 

dz 

invariant, Inv(C), can be associated with the flux through the contour C on the plane (x, y): 

Iyiy{C) = lnY^^^y){C) = £ VC(x, y)ndr = j^^ ^ VC(a;, y)^{s)ds (IV.3) 

where: n is the unit vector normal to the curve C, dr = e^dx + Gydy on the plane {x,y); 
dv 

v(s) = — denotes the "velocity" along the trajectory; and ds stands for the differential 

ds .m— ^ 
path length. Simple transformations used in Eq.( [lV.3D are: (a) ndr = e^dy — Gydx = drx u; 

(b) V({x, y){dr x u) = {u x V({x, y))dr, where u = (0, 0, 1) is the unit vector normal to the 

plane (x, y). 
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The vector product 



A(a;,y) = z/x VC(a;,y) (IV.4) 

can be considered a non-abelian generalization of the vector potential of a solenoidal "mag- 
netic field" normal to the plane {x, y) and crossing it in the points {xi.yi) and (0:2, y2)- Thus, 
A defines the flat connection of the double punctured plane z with respect to the action of 
the group G. 

It is easy to show how the basic formulae (|IV.2| ) and (|IV.3|) transform in case of commu- 



tative group Gcomm{5'i5 fl'2} which distinguishes only the classes of homology of the contour 
C with respect to the removed points on the plane. The corresponding conformal transform 
is performed by the function (^{z) = \n{z — zi) + ln{z — Z2). This immediately gives the 
abelian connection and the Gauss linking number as a topological invariant: 



A(r) = z/ X 



i={i,2} i'' 



\2' 



where ni and n2 are the winding numbers of the path C around the points Mi and M2 of 
the plane (x, y). 

Substituting Eq.( [iV.l| ) written in the Euclidean coordinates (x, y) for Eq.( |lV.30 and using 



the Fourier transform for the 5-function, we can rewrite equation ( [IV. 1|) as follows 



where 



P(^, L) = ^ r e-"^^Piq, L)dq (IV.5) 

The function P(g, L) coincides with the Green function P(ro, r = fq, L) of the non- 
stationary Schrodinger-like equation for the free particle motion in a "magnetic field" with 
the vector potential (|IV.4| ): 

'^-P(ro, r, g, L) - f ^ V - iqA{v)]\{v^, r, g, L) = 6{L)6{r - fq) (IV. 7) 



where q plays a role of a "charge" and the magnetic field is considered transversal, i.e. 
rotA(r) = 0. 

Describe now the constructive way of getting the desired conformal transform. The 
single- valued inverse function z{Q = C~^(^) is defined in the fundamental domain of ( — the 
triangle ABC. The multivalued function (f){() is determined as follows: 

- The function (f){() coincides with z{() in the basic fundamental domain; 

- In all other domains of the covering space ( the function is analytically continued 
through the boundaries of these domains by means of fractional transformations consistent 
with the action of the group G. 
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Consider two basic contours Pi and P2 on ( being the conformal images of the contours Ci 
and C2 (fig.pTlb). The function (f){z) {z 7^ {zi, Z2, 00}) obeys the following transformations: 



Ci 

z ^ z 



ai(f){z) + hi 
Ci(j){z) + di 



where 



C2 
z z 



>2U 



a2(^{z) + 62 

C2<^{z) + d2 



(IV.8) 



ai 


M 




(0-2 &2 \ 


Cl 


di 


= 91, 


\C2 d2 



92 



(IV.9) 



are the matrices of basic substitutions of the group G{gi,g2}. 

We assume ^(2;) to be a ratio of two fundamental solutions, Ui{z)^ and, U2{z), of some 
second order differential equation with peculiar points {zi = (0,0), 2:2 = (0,1), 2:3 = (c)o)}. 
As it follows from the analytic theory of differential equations jH^, the solutions ui{z) and 
U2{z) undergo the linear transformations when the variable z moves along the contours Ci 
and C2: 



Cl : 



Ui[Z) 
U2{z) 



91 



Ui[Z) 
U2iz) 



Co 



Ui[Z) 
U2{z) 



92 



Ui[Z) 
U2{z) 



(IV.IO) 



The problem of restoring the form of differential equation knowing the monodromy ma- 
trices gi and g2 of the group C known as Riemann-Hilbert problem has an old history ||68|| . 
In our particular case we restrict ourselves with the well investigated groups r2 (the free 
group) and PSL{2,'Z) (the modular group). ( [111.16| ). Thus, we have the following second- 
order differential equations: 



.(z-1)— . 



'^f\z) + (2z-l)^u^f\z) + -u 
dz 4 



(/) 



{z) = 



for the free group and 

z{z 



l)-^«(-)| 
dz^ 



z)+i-z 







(IV.ll) 



(IV. 12) 



for the modular group. 

The function which performs the conformal mapping of the upper half-plane Im^; > 
on the fundamental domain (the curvilinear triangle ABC) of the universal covering ( now 
reads 



C(^) 



u 



(f,m). 



(IV.13) 



where Ui'^^\z) and Ui^^^'{z) are the basic solutions of ( [IV.11|) and ( |IV. 12| ) for r2 and 
PSL(2,'Z) respectively. 

As an example we give an explicit form of the complex potential A{z) for the free group 
r2. Substituting Eq.(|IV3) for Eq.(|Vl3D, we get 



A{z) 



dCjz) 
dz 



'Fi{z)F,iz) F,{z) 



2(z-r 



Fi{z) 



F2(z) 



(IV.14) 
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where ^ 



The asymptotic of ([IV.14|) is as follows 

dC{z) 



z 

dz ' \ 1 



z - 1 



z 
z 1 



(compare to the abelian case). 

B. Random walk on double punctured plane and conformal field theory 

The geometrical construction described in the previous section is evidently related to 
the conformal field theory. In the most direct way this relation could be understood as 
follows. The ordinary differential equations Eq.( |lV.lI| ) and Eq.( |lV.l^ ) can be associated 
with equations on the four-point correlation function of some (still not defined) conformal 
field theory. The question remains whether it is always possible to adjust the central charge 
c of the corresponding Virosoro algebra and the conformal dimension A of the critical theory 
to the coefficients in equations like ( IV.11| ), ( |1V.12|) . The question is positive and we show 



that on the example of the random walk on the double punctured plane with the monodromy 
of the free group. 

We restrict ourselves to the "critical" case of infinite long trajectories, i.e. we suppose 
L — > cxD. In the field-theoretic language that means the consideration of the massless free 
field theory on z. Actually, the partition function of the selfintersecting random walk on z 
written in the field representation is generated by the scalar Hamiltonian H = ^(V(f)'^+mip'^ 
where the mass m functions as the "chemical potential" conjugated to the length of the path 
(m ~ 1/L). Thus, for L ^ oo we have rric = which corresponds to the critical point in 
conformal theory . 



We introduce the conformal operator, <f{z), on the complex plane z. The dimension. A, 
of this operator is defined from the conformal correlator 

MzMz')) ~ ^ (IV.15) 

\z — z \ 

Let us suppose (p{z) to be a primary field, then the four-point correlation function 
{(p{zi)(p{z2)(p{z3)ip{z4)) satisfies the equation following from the conformal Ward identity 



|69|j70| . In form of ordinary Riemann differential equation, Eq. ( |I V. 1"5| ) on the conformal 



correlator ijj{z\zi, Z2, z^) = {(p{z)(p{zi)(p{z2)(p{z3)) with the fixed points {zi = (0,0), Z2 
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1^1, 0), = 00} reads []65| , |69 

fS dp I d 1 d A A 2A 

+ + 7 7TTT + 



[2(2A + l)dz2 zdz z-ldz z^ {z - 1)^ z{z - I) \ 

1p{z\zi,Z2,Z3) = 

Performing the substitution 

IpizlZi, Z2, Z3) = [z{z - 1)]"^^ U{z) 

we get the equation 

z{z - l)u"{z) - -(1 - 4A)(1 - 2z)u'{z) - -(2A - 8A^)u{z) = (IV.16) 
3 3 

which coincides with Eq.( |IV.l"l] ) for one single value of A 

A = (IV.17) 

The conformal properties of the stress-energy tensor, T{z), are defined by the coefficients, 
Ln, in its Laurent expansion, 

00 T 

n=— 00 



These coefficients form the Virosoro algebra |]65 



[Ln, L„i\ — (n — m)Ln+m + Y2^^^'^ ~ 1T')^n+mfl 

where the parameter, c, is the central charge of the theory. Using the relation c = '^^2A+-^^ 
established in and Eq. (|IV.17| ) we obtain 

c=-2 (IV.18) 

We find the following fact, mentioned by B. Duplantier, very intriguing. As he has 
pointed out, the value A = — | (Eq.( [lV.r7|) ) coincides with the surface exponent (i.e. with 
the conformal dimension of the two point correlator near the surface) for the dense phase of 
the 0{n = 0) lattice model (or, what is the same, for the Potts model with q = 0) describing 
statistics of the so-called "Manhattan random walks" (known also as "dense polymers" — 
see the paper [^). Recall that Potts model has been already mentioned in the Chapter 
1 in connection with construction of algebraic knot invariants. It is hard to believe that 
such coincidence is occasional and we hope that the relation between these problems will be 
elucidated in the near future. 

The conformal invariance of the random walk ||66| , |67[| together with the geometrical in- 
terpretation of the monodromy properties of the four-point conformal correlator established 
above enable us to express the following assertion: 

The critical conformal field theory characterized by the values c = —2 and A = — | gives 
the field representation for the infinitely long random walk on the double punctured complex 
plane. 
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With respect to the four-point correlation function, we could ask what happens with the 
gauge connection Aj{z) if the argument Zj of the primary field ip{zj) moves along the closed 
contour C around three punctures on the plane. From the general theory it is known that 



Aj{z) can be written as 



•4.W = fE|^ (IV.19) 

where k is the level of the corresponding representation of the Kac- Moody algebra and Ri, 
Rj are the generators of representation of the primary fields ip{zi), ip{zj) in the given group 

The holonomy operator x(C) associated with Aj{z) reads 

= Pexp {-j^Aj{z)d^ (IV.20) 

It would be interesting to compare Eq. (|IV.14|) (with one puncture at infinity) to Eq. (|IV.19| ). 
Besides we could also expect that Eq. ( |1V.2D would allow us to rewrite the holonomy operator 



( [TOOD as follows 

= exp (Cin - Cfin) 

At this point we finish the brief discussion of the field-theoretical aspects of the geomet- 
rical approach presented above. 



C. Statistics of random walks with topological constraints in the two dimensional 

lattice of obstacles 

The conformal methods can be applied to the problem of calculating the distribution 
function for random walks in regular lattices of topological obstacles on the complex plane 
w = u + iv. Let the elementary cell of the lattice be the equal-sided triangle with the side 
length c. 

Introduce the distribution function P{wo,w, L\hom) defining the probability of the fact 
that the trajectory of random walk starting at the point Wq comes after "time" L to the 
point w and all paths going from Wq to w belong to the same homotopy class with respect to 
the lattice of obstacles. Formally we can write the diffusion equation 

^A^P{w, L|hom) = g^P{w, L|hom) (IV.21) 

with initial and normalization conditions: 

P{w,L = O|hom) = 5{zo); 

y2 P{wo,w, L\hom) = exp 
ru 1 vraL 

jnom} 

The conformal methods can be used to find the asymptotic solution of Eq. (|IV.21|) when 
L ^ a. Due to the conformal invariance of the Brownian motion, the new random process 
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in the covering space will be again random but in the metric-dependent "new time". In 
particular, we are interested in the probability to find the closed path of length L to be 
unentangled in the lattice of obstacles. 

The construction of the conformal transformation C(^) (explicitly described in [^) can 
be performed in two steps — see fig.|Tl|: 

1. First, by means of auxiliary refiection w{z) we transfer the elementary cell of the 
w-plane to the upper half-plane of the Im(z) > of the double punctured plane z. The 
function w{z) is determined by the Christoffel-Schwarts integral 

w{z) = ^f-Y , (IV.22) 



where -B^,^^ is the Beta-function. The correspondence of the branching points is as 
follows: 

A{w = 0) ^ A{z = 0) 

B{w = c) ^ 13{z = 1) 

C (w = c e"*i^ —>■ C{z = 00) 

2. The construction of the universal covering ( for the double punctured complex plane 
z is realized by means of automorphic functions. If the covering space is free of obstacles, 
the corresponding conformal transform should be as follows 

-^UiC)} = o (IV.23) 

where {^(C)} is the so-called Schwartz's derivative 



It is well known in the analytic theory of differential equations |68] that the solution of 



Eq.( |IV.23D can be represented as ratio of two fundamental solutions of some second order 
differential equation with two branching points, namely, of Eq. ( PV.1"T]) . The final answer 
reads 

z{C) = k\0 = ' . , (IV.24) 

where ^2(0, C) and ^3(0, Q are the elliptic Jacobi Theta-functions. We recall their definitions 

00 



^2 (x, e^"^) = 2e'-^^ e^"«"("+^) cos(2n + l)x 

n=0 

00 

^3 e^"'^) = 1 + 2 ^ e^"'^"' cos2nx 



(IV.25) 



n=0 
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The branching points A, B, C have the images in the vertex points of zero-angled 
triangle lying in the upper half-plane of the plane C,. We have from Eq. ([lV.24| ): 

A{z = 0) -> A{C = oo) 
B{z = 1) ^ 5(C = 0) 
^(^ = 00) ^ ^(( = -1) 

The half-plane Im((^) > functions as a covering space for the plane w with the regular 
array of topological obstacles. It does not contain any branching point and consists of the 
infinite set of Riemann sheets, each of them having form of zero-angled triangle. These 
Riemann sheets correspond to the fibre bundle of w. 

The conformal approach gives us a well defined nonabelian topological invariant for the 
problem — the difference between the initial and final points of the trajectory in the covering 
space (see Section 3.1). Thus, the diffusion equation for the distribution function P(C, L) in 
the covering space C with given initial point Co yields 



" p(C,Co,L) = |^'(C)|'.^P(C,Co,^) 



dL 



(IV.26) 



where we took into account that under the conformal transform the Laplace operator is 
transformed in the following way 



A, 



dw 



A, 



and 



dC 



dw 



1 



k'(c)r 



In particular, the value P{( = (q,(o,L) gives the probability for the path of length L to be 
unentangled (i.e. to be contractible to the point) in the lattice of obstacles. 

The expression for the Jacobian |w'(C)|^ one can find using the properties of Jacobi 



Theta-functions ||72|. Write w'{() = w'{z) z'{(), where 



w'{z) 



9^ 



16/3 



n8/3 ^8/3 



and 



/)4 a4 



(we omit the arguments for compactness). 
The identity 



TT 



d , (92 



d9^{x,e'<) 



x=o 



n9o{x,e'''^)92{x,e'''^)9,{x,e'''^] 



enables us to get the final expression 



\w'{0\' = c'e\9[ (O 



8/3 



h 



0.129 



(IV.27) 
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where 



1 \ n inn(n+l)C 



sm(2n + l)x 



(IV.28) 



n=0 



Return to Eq. ( |IV.26| ) and perform the conformal transform of the upper half-plane Im^ > 
to the interior of the unit circle on the complex plane r in order to use the symmetry 
properties of the system. It is convenient to choose the following mapping of the vertices of 
the fundamental triangle ABC 



A(C = oo) - 

^(C = o) - 
C{c = -i) - 

The corresponding transform reads 

•TT T 



C'iC 



e 3 
e* 3 ) 



CM 



•22r 

e 3 



e 3 



and the Jacobian |t(7'(r)p takes the form 



\w\t)[ 



T 



8/3 



(IV.29) 



(IV.30) 



In fig.(|T2]) we plot the function g{r,ilj) = — |w'(r)p where r 



r e 



iip 




FIG. 12. Relief of the function g{r, tp) — see explanations in the text. 
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The gain of such representation becomes clear if we average the function g{r,ip) with 
respect to ip- The numerical calculations give us: 



1 /•2'^ 

lim {g{r, ip)), = lim — / g{r, ip)dip = — 
where w ~ 0.0309 (see the fig.(|l3D). 



(IV.31) 




FIG. 13. Plot of product 27r {g{r, V'))^ ^ (1 ~ ^ ) ^ function of r. 

Thus it is clear that for r rather close to 1 the diffusion is governed by the Laplacian 
on the surface of the constant negative curvature (the Lobachevskii plane). Representation 
of the Lobachevskii plane in the unit circle and in the upper half-plane (i.e. Poincare and 
Klein models) has been discussed in Section 2.4. Finally the diffusion equation (|IV.26| ) takes 
the following form: 



9iV 



(IV.32) 



where P 



is the "diffusion coefficient" in the Lobachevskii plane and N = L/a is the 



dimensionless chain length (i.e. effective number of steps). 



1 + r 



Changing the variables (r,ip) (/U, '^), where /i = In , we get the unrestricted 

1 — r 

random walk on the 3-pseudosphere (see Eq.( [111.55| )). Correspondingly the distribution 
function P(/x, N) reads 



P(/x,iV) 



e 4 



^exp(-4fe) 

A'kJ2tt{NVY -^m Vcosh^ - cosh/i 



(IV.33) 



The physical meaning of the geodesies length on 3-pseudosphere, /i, is straightforward: /j, 
is the length of so-called "primitive path" in the lattice of obstacles, i.e. length of the shortest 
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path remaining after all topologically allowed contractions of the random trajectory in the 
lattice of obstacles. Hence, /i can be considered a nonabelian topological invariant, much 
more powerful than the Gauss linking number. This invariant is not complete except one 
point /i = where it precisely classifies the trajectories belonging to the trivial homotopic 
class. Let us note that the length rj is proportional to the length of the primitive (irreducible) 
word written in terms of generators of the free group 

V. PHYSICAL APPLICATIONS. POLYMER LANGUAGE IN STATISTICS OF 

ENTANGLED CHAIN LIKE OBJECTS 

Topological constraints essentially modify the physical properties of statistical systems 
consisting of chain-like objects of completely different nature. It should be said that topolog- 
ical problems are widely investigated in connection with quantum field and string theories, 
2D-gravitation, statistics of vortexes in superconductors and world lines of anyons, quan- 
tum Hall effect, thermodynamic properties of entangled polymers etc. Modern methods 
of theoretical physics allow us to describe rather comprehensively the effects of nonabelian 
statistics on physical behavior for each particular referred system; however, in our opinion, 
the following general questions remain obscure: 

How does the changes in topological state of the system of entangled chain-like objects 
effect their physical properties? 

How can the knowledge accrued in statistical topology be applied to the construction of 
the Ginzburg-Landau-type theory of fluctuating entangled chain-like objects? 

In order to have representative and physically clear image for the system of fluctuating 
chains with the full range of nonabelian topological properties it appears quite natural to 
formulate general topological problems in terms of polymer physics. It allows us: to use a ge- 
ometrically clear image of polymer with topological constraints as a model corresponding to 
the path integral formalism in the field theory; to advance in investigation of specific physical 
properties of biological and synthetical polymer systems where the topological constraints 
play a significant role. 

For physicists the polymer objects are attractive due to many reasons. First of all, 
the adjoining of monomer units in chains essentially reduces all equilibrium and dynamic 
properties of the system under consideration. Moreover, due to that adjoining the behavior 
of polymers is determined by the space-time scales larger than for low-molecular-weight 
substances. The chain-like structure of macromolecules causes the following pecuharities 
(see, for instance, [^): the so-called "linear memory" (i.e. fixed position of each monomer 
unit along the chain); the low translational entropy (i.e. the restrictions on independent 
motion of monomer units due to the presence of bonds); large space fluctuations (i.e. just a 
single macromolecule can be regarded as a statistical system with many degrees of freedom) . 

It should be emphasized that the above mentioned "linear memory" leads to the fact that 
different parts of polymer molecules fluctuating in space can not go one through another 
without the chain rupture. For the system of non-phantom closed chains this means that 
only those chain conformations are available which can be transformed continuously into 
one another, which inevitably give rise to the problem of knot entropy determination (see 
Section 2 for details). 
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A. Polymer chain in 3D-array of obstacles 



The 3D-model "polymer chain in an array of obstacles" (PCAO) can be defined as follows 
( Suppose a polymer chain of length L = Na is placed between the edges of the 

simple cubic lattice with the spacing c, where and a are the number of monomer units 
in the chain and the length of the unit correspondingly. We assume that the chain cannot 
cross ("pass through") any edges of the lattice. 

The PCAO-model can be considered as the basis for a mean-field-like self-consistent 
approach to the major problem of entropy calculation of ensembles of strongly entangled 
fiuctuating chains. Namely, choose the test chain, specify its topological state and assume 
that the lattice of obstacles models the effect of entanglements with the surrounding chains 
("background"). Neglecting the fluctuations of the background and the topological con- 
straints which the test chain produces for itself, we lose information about the correlations 
between the test chain and the background. However even in the simplest case we arrive 
at some nontrivial results concerning statistics of the test chain caused by topological in- 
teractions with the background. This means that for the investigation of properties of real 
polymer systems with topological constraints it is not enough to be able to calculate the 
statistical characteristics of chains in lattices of obstacles, but it is also necessary to be able 
to adjust any speciflc physical system to the unique lattice of obstacles, which is much more 
complicated task. 

So, let us take a closed polymer chain without volume interactions (i.e. a chain with 
selflntersections) in the trivial topological state with respect to the 3D lattice of obstacles. 
It means that the chain trajectory can be continuously contracted to the point. It is clear 
that because of the obstacles, the macromolecule will adopt more compact conformation 
than the standard random walk without any topological constraints. 

It is convenient to begin with the lattice realization of the problem. In this case the 
polymer chain can be represented as a closed A^-step random walk on a cubic lattice with 
the length of elementary step a being equal to the spacing of the array of obstacles, c. The 
general case a ^ c will be considered later. 

The random walk on a 3D-cubic lattice in the presence of the regular array of topological 
constraints produced by uncrossible strings on the dual lattice is equivalent to the free random 
walk on the graph — the Cayley tree with the branching number z = 6. 



The the average space dimension -R(A^) = yR'^{N) of the closed unentangled A^-step 
random walk is ( [0): 

R ~ aN^/^ (V.l) 

The outline of the derivation of the result (|V. 1|) is as follows. First of all note that the 

Cayley tree with z branches (called latter as z-tree with z = 2D branches) plays a role 

of the universal covering and is just a visualization of the free group Too with the infinite 

number of generators. At the same time V^jjP = where Gammaz/2 is the free group 

with z generators. Writing down the recursion relations for the probability P{k, N) for the 

A^-step random walk on the 2;-tree (compare to ( |III.24|) -( [III.53| )), we can easily find the 

T • II--- 1- -1 - r 1 r - 7~./7 I Ar\ P (k , m) P (k , N — m) 
conditional limiting distribution for the function P[k,m\i\j = -7— . Recall 

z{z — 1)'^^^ 

that P{k, m\N) gives the conditional probability distribution of the fact that two sub-chains 
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Ci and C2 of lengths m and N — m have the common primitive path k under the condition 
that the composite chain Ci C2 of length is closed and unentangled in regard to the 
obstacles; 

/ N V^'^ ( k'^N \ 

P(k, m\N) ~ — e exp ^ (V.2) 

^' ' ^ \2m{N-m)J ^\ 2m{N - m) J ^ ' 

This equation enables us to get the following expressions for the mean length of the primitive 
path, {k{m)) of closed unentangled iV-link chain divided into two parts of the lengths m and 
N — m correspondingly 



N 



k=0 




{k{m)) = Y: kPik,m\N) ^ ^^/ 2m(iV^ m) ^ ^^^^^ 



The primitive path itself can be considered a random walk in a 3D space with restriction that 
any step of the primitive path should not be strictly opposite to the previous one. Therefore 
the mean-square distance in the space ((fq — r-m)^) between the ends of the primitive path 
of k{m) steps is equal to 



;ro - r„)n = ^ka' (V.4) 

where is the radius- vector of a link with the number m and the boundary conditions are: 
r^v = To = 0. The mean-square gyration radius, (^Rl) of iV-step closed unentangled random 
walk in the regular lattice of obstacles reads 



1 1 ^ 

n^m m=l g-j 

= aViV 

z-2 8 

This result should be compared to the mean-square gyration radius of the closed chain 
without any topological constraints, <^-Rg,o) ~ 'h'^^-^ 

The relation R ~ N^^'^ is reminiscent of the well-known expression for the dimension of 
randomly branched ideal macromolecule. The gyration radius of an ideal "lattice animal" 
containing N links is proportional to iV^/^. It means that both systems belong to the same 
universality class. 

Now we turn to the mean-field calculation of the critical exponent u of nonselfintersecting 
random walk in the regular lattice of obstacles |]75|. Within the framework of Flory-type 
mean-field theory the nonequilibrium free energy, F{R), of the polymer chain of size R with 
volume interactions can be written as follows 



F{R) = F,^,{R) + F,i{R) (V.6) 

where Fint(-R) is the energy of the chain self-interactions and Fei{R) is the "elastic", (i.e. 
pure entropic) contribution to the total free energy of the system. Minimizing F{R) with 
respect to R for fixed chain length, L = Na, we get the desired relation R ~ N'^ . 
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Write the interacting part of the chain free energy written in the virial expansion 

FintiR) = V (Bp" + Cp"") (V.7) 

where V ^ R''' is the volume occupied by the chain in ci-dimensional space; p = y is the 
chain density; B = and C = const > are the two- and three- body interaction 

constants respectively. In the case B > third virial coefficient contribution to Eg. ( |V . ?[ ) 
can be neglected |73 |. 

The "elastic" part of the free energy Fei{R) of an unentangled closed chain of size R and 
length Na in the lattice of obstacles can be estimated as follows 

Fei{R) = const + In P{R, m = N/2, N) = const + In J dk P{k, N) P{R, k) (V.8) 

where the distribution function P{k, m, N) is the same as in Eq.( |V.2D and P{R, k) gives the 
probability for the primitive path of length k to have the space distance between the ends 
equal to R: 



(V.9) 



Substituting Eq.( [V.8|) for Eqs.( |V.2| ) and (|V.9| ) we get the following estimate 

/ p4 \ 1/3 



Equations ( |V.7| ) and ( |V.10| ) allow us to rewrite Eq. ( |V.6| ) in the form 

nR)^B-^-[-^) (V.ii) 

Minimization of Eq. (|V.11| ) with respect to R for fixed yields 

R ^ i?3/(4+3i?)(^c)=^/(^+3D)^.. ^ ^ (y_^2) 

The upper critical dimension for that system is D^.^ = 8. For D = 3 Eq.( |V.l^ ) gives 

R^N'^/^^ (V.13) 

It is interesting to compare Eq. (|V.12|) to the critical exponent z/an of the lattice animal 
with excluded volume in the D-dimensional space, z/an = which gives z/an = | for D = 3. 
The difference in exponents signifies that the unentangled ring with volume interactions and 
the nonselfintersecting "lattice animal" belong to different universality classes (despite in 
the absence of volume interactions they belong to the same class). 



B. Collapsed phase of unknotted polymer 

In this Section we show which predictions about the fractal structure of a strongly col- 
lapsed phase of unknotted ring polymer can be made using the concept of "polymer chain 
in array of obstacles" . 
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1. "Crumpled globule" concept in statistics of strongly collapsed unknotted polymer loops 

Take closed nonselfintersecting polymer chain of length in the trivial topological state 
0. After a temperature decrease the formation of the collapsed globular structure becomes 



thermodynamically favorable [Q. Supposing that the globular state can be described in 
the virial expansion we introduce as usual two- and three-body interaction constants: B = 
< and C = const > 0. But in addition to the standard volume interactions we 
would like to take into account the non-local topological constraints which obviously have a 
repulsive character. In this connection we express our main assertion [f79[] . 

The condition to form a trivial knot in a closed polymer changes significantly all ther- 
modynamic properties of a macromolecule and leads to specific non-trivial fractal properties 
of a line representing the chain trajectory in a globule. We call such structure crumpled 
(fractal) globule. We prove this statement consistently describing the given crumpled 
structure and showing its stability. 

It is well-known that in a poor solvent there exists some critical chain length, g*, de- 
pending on the temperature and energy of volume interactions, so that chains which have 
length bigger than g* collapse. Taking long enough chain, we define these (yf*-link parts as 
new block monomer units (crumples of minimal scale). 

Consider now the part of a chain with several block monomers, i.e. crumples of smallest 
scale. This new part should again collapse in itself, i.e. should form the crumple of the 
next scale if other chain parts do not interfere with it. The chain of such new sub-blocks 



(crumples of new scale) collapses again and so on till the chain as a whole (see fig.|T4D forms 
the largest final crumple. 



The fact that the closed chain cannot intersect itself causes two types of interactions: a) volume 
interactions which vanish for infinitely thin chains and b) topological constraints which remain 
even for chain of zero thickness. 
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FIG. 14. (a)-(c) Subsequent stages of collapse; (d) Self-similar structure of crumpled globule 
segregated on all scales. 



Thus the procedure is completed when all initial links are united into one crumple of 
the largest scale. It should be noted that the line representing the chain trajectory obtained 
through the procedure described above resembles the SD-analogue of the well known self- 
similar Peano curve. 

The specific feature of the crumpled globule is in the fact that different chain parts 
are not entangled with each others, completely fill the allowed volume of space and are 
"collapsed in themselves" starting from the characteristic scale g* . It may seem that due 
to space fluctuations of the chain parts all that crumples could penetrate each others with 
the loops, destroying the self-similar scale-invariant structure described above. However 
it can be shown on the basis of PCAO-model that if the chain length in a crumple of 
an arbitrary scale exceeds then the crumples coming in contact do not mix with each 
other and remain segregated in space. Recall that N^, is the characteristic distance between 
neighboring entanglements along the chain expressed in number of segments and, as a rule, 
the values of iVe lie in the range 30 -i- 300 |j73[ . 



Since the topological state of the chain part in each crumple is fixed and coincides with 
the state of the whole chain (which is trivial) this chain part can be regarded as an unknotted 
ring. Other chain parts (other crumples) function as effective lattice of obstacles surrounding 
the "test" ring — see fig.ll5|. 
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FIG. 15. (a) Part of the closed unknotted chain surrounded by other parts of the same chain; (b) 
Unentangled ring in lattice of obstacles. The obstacles replace the effect of topological constraints 
produced by other part of the same chain. 

Using the results of the Section 5.1 (see Eq.( |V.5|) ) we conclude that any M-link ring 
subchain without volume interactions not entangled with any of obstacles has the size 
i?(o)(M) ~ aM^/^. If R^^^ is the size of an equilibrium chain part in the lattice of ob- 
stacles, the entropy loss for ring chain, S, as a function of its size, R, reaches its maximum 
for R ~ R^^^ (see Eq. |V.10| ) and the chain swelling for values of R exceeding R^'^^M) is 
entropically unfavorable. At the same time in the presence of excluded volume the following 
obvious inequality must be fulfilled R{M) ~ aM^^^, which follows from the fact that density 
of the chain in the globular phase p ~ R^/N is constant. In connection with the obvious 
relation R{M) > R^^^ (M) we conclude that the swelling of chains in crumples due to their 
mutual inter-penetration with the loops does not result in the entropy gain and, therefore, 
does not occur in the system with finite density. It means that the size of crumple on each 
scale is of order of its size in dense packing state and the crumples are mutually segregated 
in space. These questions are discussed in details in the work ^ . 

The system of densely packed globulized crumples corresponds to the chain with the 
fractal dimension Df = 3{Df = 3is realized from the minimal scale, g*, up to the whole 
globule size). The value g* is of order 



g* = Neipa^ 



(V.14) 



where p is the globule density. This estimation was obtained in fT^ using the following 
arguments: g = (pa^)"^ is the mean length of the chain part between two neighboring 
(along the chain) contacts with other parts; consequently Nf.g is the mean length of the chain 
part between topological contacts (entanglements). Of course, as to the phantom chains, 
Gaussian blobs of size g are strongly overlapped with others because pair contacts between 
monomers are screened (because of so-called 6'-conditions |7^). However for nonphantom 
chains these pair contacts are topologically essential because chain crossings are prohibited 
for any value and sign of the virial coefficient. 

The entropy loss connected with the crumpled state formation can be estimated as 
follows: 
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(V.15) 



Using Eq.( |V.15| ) the corresponding crumpled globule density, p, can be obtained in the 
mean-field approximation via minimization of its free energy. The density of the crumpled 
state is less than that of usual equilibrium state what is connected with additional topological 
repulsive-type interactions between crumples: 

^"""'P " 1 + const(lvCiV,) ^ ^^-^^^ 

where peq is the density of the Lifshits' globule. 

The direct experimental verification of the proposed self-similar fractal structure of the 
unknotted ring polymer in the collapsed phase meets some technical difficulties. One of the 
ways to justify the "crumpled globule" (CG) concept comes from its indirect manifestations 
in dynamic and static properties of different polymer systems. The following works should 
be mentioned in that context: 

1. The two-stage dynamics of collapse of the macromolecule after abrupt changing of the 
solvent quality, found in recent light scattering experiments by B. Chu and Q. Ying (Stony 
Brook) g. 

2. The notion about the crumpled structure of the collapsed ring polymer allowed to ex- 
plain MM the experiments on compatibility enhancement in mixtures of ring and linear chains 



|85| , as well as to construct the quantitative theory of a collapse of A^-isopropilacrylamide 



gel in a poor water |83 



The paper [82] where the authors claim the observation of the crumpled globule in 



numerical simulations. 



2. Knot Formation Probability 

We can also utilize the CG-concept to estimate the trivial knot formation probability for 
dense phase of the polymer chain. Let us repeat that the main part of our modern knowledge 
about knot and link statistics has been obtained with the help of numerical simulations based 
on the exploiting of the algebraic knot invariants (Alexander, as a rule). Among the most 
important results we should mention the following ones: 

- The probability of the chain self-knotting, p{N), is determined as a function of chain 
length under the random chain closure In the work |]8^ (see also the recent paper 

) the simulation procedure was extended up to chains of order A^ ~ 2000, where the 



exponential asymptotic of the type 

po ~ exp{-N/No{T)) 

has been found for trivial knot formation probability for chains in good and 6'-solvents. 
A statistical study of random knotting probability using the Vassiliev invariants has been 
undertaken in recent work . 



The knot formation probability p is investigated as a function of swelling ratio a {a < 1) 
where a = \ (Rg) / (R'^ q), \^g) mean-square gyration radius of the closed chain 
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and y-Rg^o) ~ T2-^^'^ same for unperturbed (a = 1) chain — see fig.0, where points 

correspond to the data of Ref. [|1|; dashed hne gives approximation in weak compression 
regime and sohd hne — the approximation based on the concept of crumpled globule. 



T ' I I I I I ' I I r 
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FIG. 16. Dependence of non-trivial knot formation probability, p on swelling parameter, a, in 
globular state. 



It has been shown that this probability decreases sharply when a coil contracts from 
swollen state with a > 1 to the Gaussian one with a = 1 and especially when it 

collapses to the globular state |I|J86|]. 



- It has been established that in region a > 1 the topological constraints are screened 
by volume interactions almost completely PU[ . 

- It has been shown that two unentangled chains (of the same length) even without 
volume interactions in the coil state repulse each other as impenetrable spheres with radius 



oforder ,/(i?2 o\ 



Return to fig.|l^, where the knot formation probability p is plotted as a function of 
swelling ratio, a, in the globular region (a < 1). It can be seen that in compression region, 
especially for a < 0.6 data of numerical experiment are absent. It is difficult to discriminate 
between different knots in strongly compressed regime because it is necessary to calculate 
Alexander polynomial for each generated closed contour. It takes of order 0{l^) operations 
(/ is the number of self-interactions in the projection). This value becomes as larger as 
denser the system. 

Let us present the theoretical estimations of the non-trivial knot formation probability 
p{a) in dense globular state [a < 0.6) based on the CG-concept. The trivial knot formation 
probability under random linear chain closure, can be defined by the relation: 

9(«) = |^, q{a) = l-p{a) (V.17) 
Zo{a) 
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where Z{a) is the partition function of unknotted closed chain with volume interactions for 
fixed value of swelling parameter, a, and Z^la) is that of "shadow" chain without topological 
constraints but with the same volume interactions. Both partition functions can be estimated 
within the framework of the mean field theory. To do so, let us write down the classic Flory- 
type representation for the free energy of the chain with given a (in equations below we 
suppose for the temperature T = 1): 

F{a) = -lnZ(a) = Fi,t(a) + F,i{a) (V.18) 

where 

Fei(a) = -S{a) (V.19) 

Here the contributions Fint(a) from the volume interactions to the free energies of unknotted 
and shadow chain of the same density (i.e. of the same a) are equivalent. Therefore, the only 
difference concerns the elastic part of free energy, F^i, or, in other words, the conformational 
entropy. Thus, the equation ( |V.17| ) can be represented in the form: 



q{a) = exp(-F(a) - Fo{a)) = exp(5(a) - So{a)) (V.20) 

According to Fixmann's calculations pil] the entropy of phantom chain 5*0 (a) (5*0(0;) = 
InZo(a)) in region a <1 can be written in the following form: 

So{a) ~ (V.21) 

In the weak compression region 0.6 < a < 1 the probability of nontrivial knotting, p{a), 
can be estimated from the expression of the phantom ring entropy (Eq. ([V.21|) ). The best 
fit of numerical data [|I| gives us 

p(a) = 1 - exp (0.6 < a < 1) (V.22) 

where Ai and Bi are the numerical constant. 

The nontrivial part of our problem is reduced to the estimation of the entropy of strongly 
contracted closed unknotted ring (a -C 1). Using Eqs.( |V.l^ ) and ( |V.15| ) and the definition 
of a we find 

S{a) ^ -^a~' (V.23) 

In the region of our interest {a < 0.6) the Q;~^-term can be neglected in comparison with 
a~^. Therefore, we the final probability estimate has the form: 

p{a) = 1- A2 exp (-^«~^) (« < 0.6) (V.24) 

where A2 and A^e are the numerical constants (their values are given below). 

The probabilities of the nontrivial knot formation, p(a), in weak and strong compression 
regions are shown in fig.^ by the dotted and solid lines respectively. The values of the 
constants are: Ai = 1.2, Bi = 0.25, A2 = 0.925, = 34; they are chosen by comparing 
Eqs.(|V;2|) and with numerical data of Ref. [|. 
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3. Quasi-Knot Concept in Collapsed Phase of Unknotted Polymers 



Speculations about the crumpled structure of strongly contracted closed polymer chains 
in the trivial topological state could be partially confirmed by the results of Chapters 1 
and 2. The crucial question is: why the crumples remain segregated in a weakly knotted 
topological state on all scales in course of chain fluctuations. To clarify the point we begin 
by defining the topological state of a crumple, i.e. the unclosed part of the chain. Of 
course, mathematically strict definition of a knot can be formulated for closed (or infinite) 
contours exclusively. However the everyday experience tells us that even unclosed rope can 
be knotted. Thus, it seems attractive to construct a non-rigorous notion of a quasi-knot for 
description of long linear chains with free ends. 



Such ideas were expressed first in 1973 by I.M. Lifshits and A.Yu. Grosberg ||9^ for the 
globular state of the chain. The main conjecture was rather simple: in the globular state 
the distance between the ends of the chain is of order R ~ aiV^^^, being much smaller than 
the chain contour length L ~ Na. Therefore, the topological state of closed loop, consisting 
of the chain backbone and the straight end-to-end segment, might roughly characterize the 
topological state of the chain on the whole. The composite loop should be regarded as a 
quasi-knot of the linear chain. 

The topological state of a quasi-knot can be characterized by the knot complexity, r), 
introduced in Section 3 (see Eq. (pi.l5D ). It should be noted that the quasi-knot concept 



failed for Gaussian chains where the large space fluctuations of the end-to-end distance lead 
to the indefiniteness of the quasi-topological state. 

Our model of crumpled globule can be reformulated now in terms of quasi-knots. Con- 
sider the ensemble of all closed loops of length L generated with the right measure in the 
globular phase. Let us extract from this ensemble the loops with ri{L) = and find the mean 
quasi-knot complexity, {rj{l)), of an arbitrary subpart of length / {l/L = h = const; < 
/i < 1) of the given loop. In the globular state the probability 7r(r) to find the end of the 

chain of length L in some point r inside the globule of volume is of order 7r(r) ~ — - 

R-^ 

being independent on r (this relation is valid when La ^ R^). So, for the globular phase 
we could roughly suppose that the loops in the ensemble are generated with the uniform 
distribution. Thus our system satisfies the "Brownian bridge" condition and according 
to conjecture of the Section 3 (Eq. (|III.68D we can apply the following scale-invariant estimate 



for the averaged quasi-knot complexity J(jf{l)) 



^^^fij)) ^ /1/2 ^ ;,l/2^1/2 25) 



This value should be compared to averaged complexity \l {ffi})) of the part of the same 
length / in the equilibrium globule created by an open chain of length L, i.e. without the 
Brownian bridge condition 



^/{^fiiyj ^ I = HL (V.26) 

Comparing Eqs.( |V.23D and ( |V.26| ) we conclude that any part of an unknotted chain in 
the globular state is far less knotted than the same part of an open chain in the equilibrium 
globule, which supports our mean-field consideration presented above. Let us stress that our 
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statement is thermodynamically reliable and is independent of kinetics of crumpled globule 
formation. 



VI. SOME "TIGHT" PROBLEMS OF THE PROBABILITY THEORY AND 

STATISTICAL PHYSICS 

Usually, in the conclusion it is accepted to overview the main results and imperceptibly 
prepare the reader to an idea how important the work itself is... We would not like go not 
by a usual way and to make a formal conclusion, because the summary of received results 
together with brief exposition of ideas and methods were indicated in the introduction and 
some incompleteness of account could only stimulate the fantasy of the reader. 

On the contrary, we will try to pay attention to some hidden difficulties, which we per- 
manently met on our way, as well as to formulate possible, yet unsolved problems, logically 
following from our consideration. Thus, we shall schematically designate borders of given 
research and shall allow the reader most to decide, whether a given subject deserves of 
further attention or not. 



A. Remarks and comments to the section 2 



1. The derivation of Eqs. (p.59|) -( [11.60 ) assumed the passage from model with short- 



range interactions to the mean-field-theory, in which all spins are supposed to interact 
with each others. From the topological point of view such approximation is unphysical 
and requires the additional verification. We believe, that the considered model could be 
investigated with help of conformal theories and renormalisation group technique in the 
case of "weak disorder", i.e. when exists the strong asymmetry in choice of vertex crossings 
on a lattice. 

2. As it was shown above, utilization of Jones topological invariant with necessity results 
in the study of thermodynamic properties of a Potts model. In the work |^ was mentioned, 
that Alexander polynomials are naturally connected to a partition function of a free fermion 
model and, hence, to an Ising model. Probably, the use of similar functional representation 
of Alexander polynomials in the frameworks of our disordered model would result in more 
simple equations, concerned with statistical properties of the Ising spin glasses. 

3. All results received in this work are sticked to a model, which is effectively two- 
dimensional, since we are interested in statistical properties of a planar projection of knot 
in which all space degrees of freedom are thrown away and the topological disorder is kept 
only. Thus, physically, the model corresponds to the situation of a globular polymer chain 
located in a narrow two-dimensional slit. In connection with that the following question 
is of significant interest: how the space fiuctuations of a trajectory in a three-dimensional 
space modify our consideration and, in particular, the answer ( |11.67| )? 



B. Remarks and comments to sections 3 and 4 



1. The investigation of topological properties of trajectories on multiconnected mani- 
folds (on planes with sets removed points) from the point of view of the conformal field 
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theory assumes a construction of topological invariants on the basis of monodromy prop- 
erties of correlation functions of appropriate conformal theories. In connection with that 
there is a question concerning the possibility of construction of conformal theory with the 
monodromies of the locally-free group considered in work. 

2. Without any doubts the question about the relation between topological invariants 
design and spectral properties of dynamic systems on hyperbolic manifolds is of extreme 
importance. The nature of mentioned connection consists in prospective dependence between 
knot invariants (in the simplest case, Alexander polynomials), recorded in the terms of a 
trace of products of elements of some hyperbolic group (see expression ( [111.14] )) and trace 
formulae for some dynamic system on the same group. 

3. Comparing distribution function of primitive paths fi ( |1II.58|) with the distribution 
function of a knot complexity of t] ( [111.15| )), we can conclude that both these invariants 
have the same physical sense: a random walk in a covering space, constructed for lattice of 
obstacles, is equivalent from the topological point of view to a random walk on a Cayley 
tree. Thus, the knot complexity is proportional to a length of the primitive (irreducible) 
word, written in terms of group generators, i.e. is proportional to a geodesic length on 
some surface of constant negative curvature. We believe, that the detailed study of this 
interrelation will appear rather useful for utilization of algebraic invariants in the problems 
concerning statistics of ensembles of fluctuating molecules with a fixed topological state of 
each separate polymer chain. 

3. Questions, considered in Section 3 admit an interpretation in spirit of spin glass 
theories, discussed at length of the Section 2. Let us assume, that there is a closed trajectory 
of length L, which we randomly drop on a plane with regular set of removed points. Let one 
point of a trajectory is fixed. The following question appears: what is a probability to find a 
random trajectory in a given topological state with respect to the set of removed points? The 
topological state of a trajectory is a typical example of the quenched disorder. According to 
the general concept, in order to find an appropriate distribution function (statistical sum), 
it is necessary to average the moments of topological invariant over a Gaussian distribution 
(i.e. with the measure of trajectories on a plane). The same assumptions are permitted 
us to assume, that the function gir^ip) (Jacobian of conformal transformation) — see Fig.p!^ 
has a sense of an ultrametric "potential", in which the random walk takes place and where 
each valley corresponds to some given topological state of the path. The closer r is to 1, the 
higher are the barriers between neighboring valleys. Thus, all reasonably long {La 3> c^) 
random trajectories in such potential will become "localized" in some strongly entangled 
state, in the sense that the probability of spontaneous disentanglement of a trajectory of 

length La is of order of exp ( —const — ). Probably this analogy could be useful in a usual 



c2 

theory of spin glasses because of the presence of explicit expression for the ultrametric Parisi 
phase space ( [|^) in terms of a double-periodic analytic functions. 



C. Remarks and comments to section 5 



1. We would like to express the conjecture (see also 0]) concerning the possibility of 
reformulation of some topological problems for strongly collapsed chains (see Section 5.4) 
in terms of integration over the set of trajectories with fixed fractal dimension but without 
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any topological constraints. 

We have argued that in ensemble of strongly contracted unknotted chains (paths) most 
of them have the fractal dimension Df = 3; 

We believe that almost all paths in the ensemble of lines with fractal dimension Df = 3 
are topologically isomorphic to simple enough (i.e. close to the trivial ) knot. 

Let us remind, that the problem of the calculation of the partition function for closed 
polymer chain with topological constraints can be written as an integral over the set Q of 
closed paths with fixed value of topological invariant (see Chapter 1): 

Z = J D^{r}e-^ = /•••/ - Jq], (VI.27) 

n 

where D^{r} means integration with the usual Wiener measure and 6[I — Iq] cuts the paths 
with fixed value of topological invariant (Jq corresponding to the trivial knots). 

If our conjecture is true, then the integration over Q in Eg. ( [VI. 27]) for the chains in 
the globular phase (i.e. when La 3> R"^) can be replaced by the integration over all paths 
without any topological constraints, but with special new measure, Df{r}: 

Z = J . . . J Df{r}e-'' (VI.28) 

The usual Wiener measure Duj{r} is concentrated on trajectories with the fractal dimen- 
sion Df = 2. Instead of that, the measure Df{r} with the fractal dimension Df = 3 for 
description of statistics of unknotted rings should be used. 

2. We believe, that the distribution of knot complexity found for some model systems 
can serve as a starting point in construction of a mean-field Ginsburg-Landau-type theory 
of fluctuating polymer chains with a fixed topology. From a physical point of view it seems 
to be important to rise the mean-field theory which takes into account the influence of 
topological restrictions on phase transitions in bunches of entangled directed polymers. 

3. Let us note, that despite a number of experimental works, indirectly testifying for 
the existence of a fractal globule (see section 5 and references), the direct observation of 
this structure in real experiments is connected to significant technical difficulties and is so 
far not carried out. We believe, that the organization of an experiment on determination 
of a microstructure of an entangled ring molecule in a globular phase could introduce final 
clarity to a question on a crumpled globule existence. 
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